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PREFACE 


This report covers work performed by the Atmospheric Radiation 
Group of the Institute of Geophysics and Planetary Physics, UCLA 
from July 1970 through June 1972 under contract NAS 1 -10107 » spon- 
sored by NASA's AAFE program. The work was carried out under the 
direction of the principal investigator, Professor Z. Sekera. 

In this preface, the editor would like to review the goals and 
progress of the "program, and problems encountered in carrying out 
these goals. Some of these problems are discussed in detail in the 
report proper, but others do not conveniently fit into the organi- 
zation of the report. 

The major objectives of the program have been to study the 
feasibility of measuring the polarization parameters of the scat- 
tered sunlight emerging from the earth's atmosphere from a satellite 
platform and to study how such measurements may be used in the study 
of atmospheric turbidity ("Atmospheric turbidity," aerosol," and 
"suspended particulates" are essent ial 1 y synonymous terms and are 
used interchangeably throughout this report.) leading naturally to 
a division into an experimental and a theoretical program. 

It is the determination of atmospheric turbidity parameters 
from the polarization measurements which constitutes the inversion 
problem, the major concern of the theoretical program. Professor 
Sekera has prepared a highly detailed discursive introduction 
(Chapter 1) to this complex problem in which he summarizes the 
theoretical developments leading to the initiation of the current 
project. 

Essential points are reiterated, risking redundancy to dispel 
obscurity. Step-by-step analytical derivations are omitted but 
complete references to material already published in the professional 
literature are given. A number of papers resulting from work done 
under this contract will also be published,' and the results and 
significance of these papers are also discussed in this report in 
Chapters 1 and 2. No justification for the interest in particulates 
in the atmosphere is given since abundant articles on the topic 
have appeared in popular, professional and trade publications in 
recent years. 

The majority of funds were devoted to the experimental program 
which involved the construction, testing, and use in flight exper- 
iments of an engineering model of a photopol ar imeter (PPM) suitable 
in principle for satellite application. The instrument was con- 
structed under sub-contract by the General Electric Space Center, 
Valley Forge, Pennsylvania in accordance with the breadboard design 
developed under an earlier sub-contract . Description of the PPM 
and detailed test reports are presented in the Appendix. Airborne 
experiments performed on board the Convair 990 (CV-990) stationed at 
NASA/Ames Research Center, Moffett Field, Ca. are reviewed in Chap- 
ter 3, together with additional discussion of the instrumentation. 

i i i 



Even though baloon exper iments. would have more adequately simulated 
satellite measurements* due to the altitude attainable, it was 
desirable to have adequate control over, the target area selection 
which balloons can seldom assure. Since the aircraft measurements 
were not at the top of the atmosphere, it was necessary to make 
calculations applicable to intermediate levels of the atmosphere. 
Specific calculation relevant to the CV-990 experiments are presen- 
ted with the data in Chapter 3- 

Two major problems in the experimental program were a three- 
month delay in delivery of the PPM and delays in obtaining suitable 
flight dates on the CV-990. The three months would have been in- 
valuable in assessing the behavior of the instrument and correcting 
the defects. The delays in flight time necessitated by commitment 
of the CV-990 to the Barium ion cloud (B1C) experiment postponed 
the UCLA experiment into a period which was undesirable in terms of 
obtaining both suitable weather and high aerosol concent rat i ons . 

An additional month was lost due to an aircraft engine failure (and 
to frequent front passages in winter season). It was a matter of 
luck that the clouds cleared over the Mexican desert on the last 
days available for flight. 

Theoretical progress was also delayed by the difficulty of 
finding someone at a post-doctoral level qualified to work on the 
theoretical problems. It was six months before Professor Sekera 
obtained the services of Dr. Jacob Kuriyan, who, within a short 
time, was making important contributions to the fundamental theory. 
(These contributions are discussed in Chapters 1 and 2.) At the 
same time. Dr. R. E. Bradbury took over management of the experi- 
mental program from Dr. N. B. Manson who left to accept an academic 
post in Great Britain. 

Other personnel and their major functions were: Mr. R. D. Kopa 

design review and inspection of the PPM engineering aspects and 
design and installation of the airborne scan platform; Mr. D. 
Friedman, interface electronics design and construction; Mr. J. 
Sherman, computer programming. The program is also indebted to 
Dr. L. L.. Stowe who often acted in lieu of and as liason with 
Dr. Sekera though he was not directly supported under the contract. 
The contract provided partial support of Dr. T. Takashima's thesis 
work and the resulting computation of the internal field. 

Though not directly supported by contract funds, the continu- 
ing work of Dr. N. C. R. Rao of the Meteorology Department (with 
the computational assistance of Dr. Takashima) is essential to the 
program since it demonstrates the usefulness of polarization 
measurements from the top of the atmosphere, and therefore is dis- 
cussed in this report. 



'The contract called for the following five items to be dis- 
cussed in the report: 

1) A description of the analytical and experimental procedures 
for incorporating the effect of terrestrial surface reflection on 
the polarization of reflected sunlight. 

2) A description of the analytical procedures for incorpor- 
ating the effect of radiative transfer in a turbid atmosphere on 
the polarization of scattered sunlight emerging from the top of 
the earth's atmosphere. 

3) An analysis of the performance of the engineering model 
together with a description of system operation. 

4) A comparison of the analytical prediction of atmospheric 
aerosol profiles based on engineering model measurements to direct 
"in situ" measurements of aerosol size and vertical distribution. 

5) A determination of the feasibility of conducting this 
experiment from an unmanned orbiting spacecraft platform. 

At the initiation of the contract, the principal investigator 
pointed out that it would be impossible to complete the theoretical 
work within the time and money allocations of the contract, even if 
the problems mentioned above had not arisen. Indeed, if the five 
items above are not discussed in the depth or degree desired by 
NASA, sufficient support should be provided to continue the research 
during the period preceeding preparation and development for actual 
satellite f 1 i ghts . 

Item 4 provides an excellent example. This item is totally 
satisfied by the discussion in Chapter 3, where comparison of the 
data with the theoretical curves for a Rayleigh atmosphere confirm 
the insignificant aerosol concentrations as measured "in situ." 

UCLA had predicted that it would be difficult to find significant 
turbidity and suitable cloud conditions during the winter, but des- 
pite NASA's promise that there would be no difficulty in obtaining 
fly.ing time, a variety of conditions forced the flights to be post- 
poned until December, 1971* 

I terns 1 and 2 are discussed in Chapter 1 and 2 to the degree 
permitted by the theoretical progress. 

I tern 3 is covered in detail in the appendix and in portions 
of Chapter 3* The result of the analysis of the performance of the 
photopol ar imeter leaves no doubt that it is possible to construct 
satellite instrumentation capable of making high quality polariza- 
tion measurements . It is also pointed out in Chapters 1 and 2 that 
several methods may be employed to determine the total atmospheric 
turbidity from such measurements. Particular attention is called 
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to Section 2.4. 

This confirms the feasibility of conducting such experiments 
from an unmanned space platform, (though they may as well be con- 
ducted from a manned platform such as the space shuttle) thereby 
answering Item 5. Additional comments on Item 5 are made ,in Chapter 
4. Any further detailed engineering analysis with respect to orbits, 
tracking, etc. should be contracted to a potential industrial 
satellite fabricator in preparation for carrying the experiment in 
space. 

Simultaneously support for theoretical work at UCLA on a 
continuing basis should be initiated to construct a library of 
suitable models, investigate how much more information may be 
obtained from the measurements and devise more tractable computer 
programming methods to make extracting information from the meas- 
urements speedier and less costly. 
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DEVELOPMENT OF AN EXPERIMENT FOR VISIBLE RADIATION 
POLARIZATION MEASUREMENTS FROM A SATELLITE 


By Zdenek Sekera and Richard E. Bradbury 
Institute of Geophysics and Planetary Physics 
University of California, Los Angeles 


INTRODUCTION 


The Inversion Problem in Turbid Atmospheres . - Chandrasekhar 
(ref. 1) has shown that, in a plane parallel atmosphere composed * 
of Ray 1 e i gh-type scatterers, the radiative transfer problem includ- 
ing multiple scattering is soluble. The intensity and polarization 
parameters of the scattered sunlight have been numerically evaluated 
(ref. 2, 3, 4) and measured values of skylight polarization under 
very clear atmospheric conditions agree with the theoretical values, 
(ref. 5, 6, 7, 8). 

Turbid atmospheres contain suspended particles which, because 
of their size and form, do not obey the Rayleigh law of scattering. 
This results in skylight polarizations which differ from those for 
pure Rayleigh atmospheres, thereby opening a new way to study at- 
mospheric turbidity. If the computation of the effect of the par- 
ticulate matter were as easy as for a pure molecular atmosphere, 
the determination of atmospheric purity and other related quanti- 
ties would be a straightforward procedure. Unfortunately, this is 
not the case; the law of scattering of light in a polydisperse 
medium is much more complicated, and the determination of particulate 
characteristics from optical measurements requires highly complex 
and sophisticated procedures. 

The main reason for this difficulty is that the characteristics 
of scattering by particulates depends on many more parameters (e.g., 
particle number-size distribution, spatial distribution, etc.) than 
can be measured. Consequently, it is impossible to use a d i rect 
inversion method, i.e., direct computation with the observed quan- 
tities alone as input data, such as is used in temperature deter- 
minations from satellite measurements. Since direct measurements 
of all the different parameters involved requires sophisticated and 
tedious techniques yielding results which may be applicable only in 
a limited region, it is necessary to resort to calculations based 
on a pre-selected set of data, a so-called "model." 

The selection of which parameters are of interest is obviously 
determined by the purpose for studying the atmospheric turbidity, 
and the larger the scale of interest, the more restricted is the 
number of parameters. Fortunately, the physical processes in the 
atmosphere considerably restrict the size ranges of the particles 
and disperse products of local origin over a relatively small, 
isolated area. The general characteristics of particles normally 
responsible for turbidity over larger areas and throughout all 



levels of the atmosphere are more or less known to the atmospheric 
physicist and he has at his disposal a rather restricted number of 
models. if the scattering behavior and therefore the polarization 
of light emerging from the atmosphere for his set or "library" of 
models has. been calculated, he can determine the atmospheric tur- 
bidity by comparing the pol ar izat i on, measurements with his library 
and selecting from this the model that matches the data. This 
process will by no means be unique. The confidence in the uniqueness 
may be increased by utilizing measurements at a number of different 
wavelengths and requiring that the model calculation agree with the 
measurements at all wavelengths. 

Since there are usually more parameters of atmospheric turbid- 
ity than equations from which they are to be determined using the 
polarization measurements, it is obvious that there is no other 
method available than that of the "library technique." No direct 
inversion method exists for the general case and most likely will 
never exist, the theory of radiative transfer in a general turbid 
atmosphere being untractable and too difficult for even the best 
existing computers. 

The reasons for these difficulties are well understood and 
are due to the interaction of the incident radiation with the 
scatterer. If the interacting matter is an electromagnet ic dipole, 
the resulting field of scattered radiation is quite simple (Rayleigh 
scattering). Complications appear and are augmented when adjacent 
dipoles are grouped into multipcles of various orientations to form 
particles of different sizes and shapes. The differing positions 
of these individual centers produce complex interference patterns 
becoming more bizarre with increasing particle size. 

It is well known that the analytical theory of scattering by 
large particles is almost untractable and attempts have been made 
to solve only a few special cases (ref. 9 ). Fortunately-, the ef- 
fect of the shape of the particle decreases with its size and the 
cases typically encountered in atmospheric studies can be approxi- 
mated well by treating the particles as small dielectric spheres 
and hence almost all mathematical analyses have treated aerosol 
scattering as Mie scattering. 

The construction of libraries for entirely general cases is 
almost impossible; not only the computation of entries but the 
selection of entries corresponding to measured values, represents 
enormous effort even for the most modern facilities. (The problem 
might be likened to thumbing through a library card catalog arranged 
by author to find a call number matching one we have.) It is imper- 
ative to introduce further approximations, inevitably restricting 
the number of parameters or character of information which may be 
deduced from the polarization measurements of the radiation emerging 
from the earth's atmosphere. ' 
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In the following we will list some such approximations and 
the type of information which may be derived from them. We can 
begin with the case' of a pure molecular (Rayleigh) atmosphere, for 
which (spectral) optical thicknesses are known or can be calculated. 
The corresponding polarization quantities can be found in' existing 
tables or calculated by relatively simple and straightforward means. 
With a library consisting of such tables for specific wavelengths, 
(nadir) angle of observation and angle of solar illumination, we 
can compare the measured value of the degree of polarization with 
the computed value for the molecular atmosphere. Any difference 
indicates the presence of some kind of turbidity. By simple inter- 
polation, we can determine the number of additional molecules (or 
the equivalent optical thickness) needed to produce the same 
polarization. The ratio of the two optical thicknesses, r observed/ 
T"theoret i ca 1 may serve as a rather primitive measure of the atmos- 
pheric turbidity. This is, however, only a qua 1 i tat i ve measure 
answering the question "Is the atmosphere turbid, or not?" By more 
sophisticated manipulations with satellite polarization data, one 
can instead of a library technique, use a d i rect i nvers ion method 
which yields either four independent values of the optical thickness 
if the earth's surface is fully absorbing (zero albedo) or two 
independent values of the optical thickness together with two inde- 
pendent estimates of the albedo for the case of Lambertian 
(isotropic, unpolarized) reflection (ref. 10). (The problem of the 
effect of ground reflection will be discussed in greater detail in 
a later sect ion. ) 

If the direction or angle of observation is preselected, then 
the inversion i s, s impl i f i ed. The optical thickness can be read 
from a simple diagram for any zenith angle of the, sun. In Fig. 

(1-1) such a diagram is shown for the case when the polarimeter is 
directed toward the nadir. The methods discussed above are based 
on the simplifying assumption that the particulates in the atmos- 
phere can be treated as molecules (or submicron particles) and their 
effect can be noticed only. in the increase of the optical thickness 
of a pure molecular atmosphere. In this way, one may arrive at an 
"effective" optical thickness of a turbid atmosphere and although 
this measure is almost purely qualitative, it is surprising in how 
many problems it may be usefully applied. 

For any more quantitative assessment of atmospheric turbidity 
one has to consider the characteristics and behavior of scattering 
by particulates that differs from molecular scattering. 

The Forward Peak in Aerosol Scattering . - While in Rayleigh 
Scattering the polarization parameters (such as those introduced by 
Stokes, discussed in a later section) depend only on the scattering 
angle (between incident and scattered ray), they depend critically 
on the size in large particle scattering. Therefore, in the models 
of turbid atmosphere, size distribution of , part i cul ates in a unit 
volume, i.e., variation of the number-density of particles of a 
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given size, enters as a very important factor. 

Size distribution of particulates in the atmosphere obviously 
depends on the origin of these particulates, and there is practic- 
ally no limit on the number of different types that could be included 
in models of our library. Fortunately, for the natural aerosol, i.e., 
part i cul ates or i g i nati ng in the atmosphere, the physical processes in 
the atmosphere are leading to a very reasonable guess of the form 
of the particle distribution, especially in haze or dust. The coa- 
lescence of the minute water or ice coated condensation nuclei is 
responsible for the increase of the number of submicron particles, 
while the sedimentation of micron and larger particles produces a 
sharp decline of the number-density of large particles. The result- 
ing curve will then be bell -shaped with the maximum around ly. 

Such a curve can be approximated by a three-parameter one of the 
form Kl ( a ) = Noa a e"ba Y where N(a) is the number of spherical particles 
of radius a in a unit volume. The size distribution can then be 
characterized by only three parameters cj, b and Y . 

The scattering characteristics of polydispersions with such 
distribution for a few selected sets of these parameters have been 
thoroughly discussed by Deirmendjian (ref. 11) who considered Mie 
particles with the above size distribution and by integrating over 
the size derived the polarization parameters of a unit volume of 
such polydispersions. In the dependence on the scattering angie 
(all these parameters are illustrated in several instructive pictures 
in his outstanding book of tables), one can notice a pronounced 
maximum in the forward direction, usually denoted as the "forward 
peak." The origin of this peak is the addition of the diffraction 
peak of each single particle in the considered volume. From this 
fact, all the characteristics of a forward peak follow, such as the 
larger the particles in the volume, the more pronounced the peak. 
Consequently, for a rain model the intensity of scattering is 103 
larger in the forward maximum than in a sideward direction. The 
forward peak is without any doubt the most important feature of the 
scattering by such polydispersions. This becomes even more obvious 
if the intensity or other parameters are plotted in a 1 i near scale, 
as was done in Fig. (1-2) by Takashima (ref. 12) for two of 
De i rmend j i an 1 s models, and not in a logarithmic scale as is usually 
done. This forward peak is also the origin of the difficulties 
occurring in the evaluation of multiple scattering in a turbid 
atmosphere. If the elements of the phase matrix, defining the re- 
lation between polarization elements of the incident and scattered 
radiation, are developed in a Legendre series of the scattering 
angle, it is necessary to include a considerable number of the terms 
in this series in order to accurately reproduce the forward peak. 

The solution of radiative transfer equations reduced then to the same 
number of separate integral equations as the number of the terms 
in the Legendre series for the elements of the phase matrix. In 
this complicated procedure, truncating errors can easily shoot up 
to a magnitude rivaling that of the resulting value. Analogous 



difficulties occur in the so-called "Monte Carlo" method. Here 
the intensities change their meaning into probab i 1 i t i es of photons, 
and in the case of a strong forward peak (of the order 10" 5 or higher) 
it is necessary to perform the computation with an abnormally large 
number of photons, if the statistics outside the peak are to give 
rel iable resul ts. 

To surmont the difficulties of radiative transfer in a turbid 
atmosphere, it is quite appropriate to look for some simplification 
or suitable approximation. The fact that the forward peak is so 
pronounced and differs clearly from the scattering at other angles, 
leads naturally to the "method of separation of the forward peak." 
This method was applied fi rst in the atmospheric problems by Romanova 
(ref. 13) who replaced the forward peak by a delta function. That, 
however, differs quite substantially from the actual intensity 
distribution in the forward region and since this peak contains the 
largest intensity, the approximation should be the best possible. 
Weinman (ref. 14) used a Gaussian distribution for a scattering 
function in a rain cloud, and gained surprisingly good agreement. 

Had he used another normalization constant, however, he undoubtedly 
would have gotten even better agreement. The forward peak of the 
dispersion or of a single particle can be approximated by innumer- 
able ways more or less successfully. It is quite clear, however, 
that the best approximation will be obtained from approximations that 
are based on physical law. This has been fully demonstrated by a 
few unsuccessful attempts by students at UCLA. The physical back- 
ground for successful approximations has been given by Saxon (ref. 

15) who from Maxwel 1 1 s equat ion derived an expression for scattered 
field as a combination of volume and surface integrals. These in- 
tegrals were then evaluated by the K.W.B. method and Deirmendjian 
(ref. 16) has shown that the values of the intensity in the forward 
direction agree very well with the Mie values computed by Penndorf 
(ref. 17) (1955)* With a si ight numerical correction (including X = 
2rra/x) a complete agreement with Mie theory could be found in the 
entire range 3 £ X 100. 

This extraordinary agreement of the theory with Mie-calculation 
motivated us to continue the previous research and start the study 
of analytical approximation of forward peak with scattering of 
single particle approximated by the diffraction law well known to 
physicists before or derived recently by Saxon. 
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Fig. (1-2) -- Normalized Scattering Intensities of the 
Two Aerosol Models (after Dei rmendj ian) . 



2.0 THEORETICAL PROGRAM 


2 . 1 The Forward Peak Approximation and the Effective Size 
Distribution for Aerosol Scattering . - It was decided that obtaining 
an analytical approximation of the forward peak was the fundamental 
task of the theoret ical part of' this . project and the first, there- 
fore, to pursue. Following De i rmend j i an ' s ( ref . 1 T) notation, the 
elements of the phase matrix Pj (i = 1, 2, 3) all indicate the 
forward peak, i.e., maximum at‘ 9= 0, wi th, decrease toward the lar- 
ger value of scattering angle 6 . As mentioned before, this feature 
is due to the addition of the diffraction peak of each single particle 
in the volume in question. Hence the analytical expression for Pj 
for small 9 ( < 10°) is the integral 

, . ’ . • ■ • .■ v - j ; - 

P:(9) = Td-O ,a)N(a)da (1-1, 2, 3) (1) • 

1 0 1 

where Dj is the diffraction peak of single sphere of radius a, and 
N(a) - the weight function - the previously mentioned number-density, 
(number of particles of radius a in the volume) given by ..the so- 
called polygamma function, indicated before, with three independent 
parameters a, b, y . Among these parameters, a and Y determine the 
size of the distribution, and b is related to the position of the. 
maximum (1/b = (aA )a max ). 

The project director has shown that by simple change of variable 
the integral P;(9) can be transformed into a Laplace integral. In 
this way.,, the existing tables of Laplace transformations, as well as 
the calculus associated with these. transformations may be used for 
evaluation of the forward peak integral for large class of functions 
Dj and of the parameters o. , b, y . 1 

A meaningful approximation in Eq. (2-1) will be obtained if 
the approximation used for Dj(9) is valid in the region where the 
given size distribution N(.a) is a maximum. • Th i s. pr i nc i p) e,. stated 
for the first time by Kuriyan* contributed greatly to the 
success of our approximation of the forward, peak. He considered 
first the Born approximation for Dj(9,a) and found that this repre- 
sented a good approximation if 2 the wei ght ^funct ion N(a) was vanish- 
ingly small for all a > X/2rr(m 4 -l), and (rr-1 )<1. Almost all models 
that Deirmendjian labels as "haze" satisfy these conditions. On the 
other hand, the Fraunhofer approximation for Dj(9,a) is. valid in 
the range Z 73 .;( rrr -. 1 ) > X(m2-1) <c 1 . Thi.s would then appl y more for 
the ; -region of raindrops. ,.T ? he idea of different validity regions 
for the approximation of diffraction by a single particle has not 
been realized and applied before, and the success of the approxima- 
tion discussed here as well- as the explanation of a few failures by 
previous authors, should be attributed to Kuriyan's ability. 

There are several ways to evaluate the integral in (2-1). 

They will be given elsewhere. If the expression for D; 

in the- Born approximation is substituted in (2-1) in the i ntegrand, 
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the square of the Bessel, function J-, 2 is developed into a series 
and integrated term by term, one gets easilty 

- , • •• V (-) n (2»n) 

Pi (0) = P 1 (0)/P 1 (0 = ,0) = l. n! ( 4, n ) ( 5/2, n ) 

, • • n=0 

J r(e + 2n/y . ) 2 k ( b/2 ) "si n0/2 2n (2) 

. ? . t( e.) 


where e (a + 7)/ and, as usual, k = 2tt/\ (k = wavenumber of 
the monochromatic radiation) (a,n)s t( a+n)/T( a ), t ( ), the Gamma 
funct ion. 

For an integer y, P\(0) reduces to the general hypergeometr ical 
function U F V *> that then for y = 1 assumes the form ^F 2 

P,(9) = 3 F 2 ( 2, zll, ( e+1 )/2 ; 4, 5/2; -z) (3) 

•with 

*/ z = 4kb ^ sin0/2., (4) 


For integer a> Pi(0) in (2-3) reduces even further either to finite 
series or a closed expression. Hence, for a =2, 



For the other matrix elements we get for the Born approximation 

P 2 (0) = cos 2 0 P" | ( 0 ) i P 3 (0) = cos ( © ) P", ( 0), P 4 (0) = 0. 

'• - ( 6 ) 


* u F v (p r p 2’ •••’ p u ; q r q 2 ,,,, ’ q v ;z ) = ^ ( p 1 > n)...(p ll ,n) 

n-0 n! 
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In the expression for P i ( 9 ) there is only one variable appearing, 
namely z - 1 6k^b “ 2 sin-9/2. Hence all distributions for a given 
pair of parameters and Y have the same form of the _forward peak 
or belong to the same curve. In Fig’. (T- 1 ) the curve P|(0) as a 
function of z is compared with the values of z computed for various 
scattering angles from De i rmend j i an ' s tables (22-23). The agree- 
ment is very good, even in a region where the intensity has dropped 
to 30% of the forward intensity ( P ] ( 9) = 0.4). The tables of 
Deirmendjian where x - ka violates the Born criterion have been de- 
liberately omitted. 

The evaluation of the integral (2-1) for Y ^ 1 is very elab- 
orate, although it can be done for any a and y • It is interesting 
to mention that for Y = 2 the expressions in the Fraunhofer approx- 
imation contain as a factor exp '[-H ka^s i n 9)S ( ao = b“ 2 ), Weinman's 
(ref. 14) empirical solution. That would then explain such a good 
agreement with De i rmend j i an 1 s data. 

In the parametrization used by Deirmendjian for his poly- 
dispersions, namely Y = 1, 2, the theoretical expressions for 

the forward peak i ntens i ty Tj ( 0) are simple for' y = 1» but rather 
complicated for Y = ?, leading to infinite series, the convergence 
of which still should be studied and established. The closed 
analytic expression can be applied to the haze H ( Y = 1 , a = 2) 
but not for the haze of the type L and M, for which Y = Kuriyan 
has correctly pointed out that, since y enters in the exponent of 
the weight function, exp[-ba Y ], it may be possible by retaining y 
= 1 to reproduce the basic features of the size distribution (within 
the limits of experimental error) that Deirmendjian obtains with y 
= i by simply changing b. This is quite justifiable because the 
experimental data are such that the curves for Y = i and y = 1 are 
barely distinguishable (especially in the critical region, i.e., near 
maximum particle concentration). Within this educated guess, Kuriyan 
has retained the values of the parameters a = 2 and y =* 1 for all 
types of haze particle distributions and then has varied the para- 
meter b to reproduce De i rmendj ian ' s results for haze types H, L, 
and M. 

To illustrate such procedure for haze M, Kuriyan first compu- 

ted from tables 12-16 (ref. 11) one value of P i ( 6 ) , say pp for fixed 
X and 0 (0 < 0 < 10°). From the analytic formula (2-3) he found 
a value z, say z\, that corresponds to p,. Since X and 0 are fixed 
and known, he determined the new value of b, say b^. This process 
had been repeated for a few points and the average value of bw 
obtained. With this val_ue b^ for all those z values from 0 <0 
< 20° from tables 12-16 P|(o) was computed and compared with original 
De i rmend j i an ' s values. Fig. (2-2) demonstrates such comparison and 
shows again that all these calculated points (from Mie theory) fall 
on the curve, showing that the analytical formula (2-3) for a H dis- 
tribution (with a new b^) reproduces in the forward direction 
De i rmendj i an ' s exact phase matrix element for haze model M (with a 
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different b). A similar analysis has been performed on the haze 
type L. The new value b[_ that leads to the same point on the curve 
for haze L when computed from (2-3) with a = 2, _y_ = 1, has been 
computed as before, and the resulting values of P](o) plotted in 
Fig. (2-3) in a curve that again represents a good approximation 
of De i rmend j i an 1 s computations (ref. 11), (Tables 19-26). 

Summarizing previous discussions, we see that significant 
scattering properties of these different types of hazes can be 
described with good accuracy by using only one type, leading to a 
very simple analytical expression: ^( 9 ) in ( 2 - 3 ) for a. = 2 and 

y = 1 . It Is quite conceivable that an assembly of different 
objects could give rise to the same scattered field. The unique 
size or shape of the constituent haze particles, however, cannot 
be concluded from these results. It has only been established 
that, in the forward direction, the scattering can be related to 
one fundamental model characterized by the weight function (number 
density) N( a )~ a^e~° a . In conclusion, it should be mentioned that 
the formulae presented imply that in the forward direction the 
normalized phase matrix elements are only functions of one variable, 
2k sinB/2. This is, therefore, a necessary condition for the 
validity of an approximation. 

This observation will help decide if the data points for 

P](B) can be approximated by oyr formulae. For instance, if r]( 0 ) 
is plotted against 2 k sinfl/ 2 , and if the points for various k and 0 
cluster around one curve, then the formula may be applicable. If, 
on the other hand, the points are distributed randomly, then our 
formula is not appl i cabl e . To illustrate this method we used 
De i rmend i i an 1 s haze M tables 1-15 where the wavelength varies be- 
tween 0.45 u, and 16. 6u and the real part of the refractive index 
varies in the vicinity of 1.3* The scattering angle 0 is restric- 
ted to be in the range 0 <c B £ 30 ° and the plot is made (see 
Fig. (2-4) of Pj(p) against 2k sinfi/2. The points, evidently, lus- 
ter around one curve and this reassures us that our approximation 
is a consistent one. 
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HAZE M (tables 1-15) 0.45/x<X< 16,6^ 


Fig. ,(2-4) -- Normalized Forward Scatter ing "as a 
Function of 0. 
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2. 2 Separation of the Forward Peak in the Radiative Transfer 
Equat i ons . - As discussed in the previous section, it is possible 
by using the Born approximation -for the diffraction by a single 
particle, and then integrating over the size range, to achieve not 
only a very good approximation wi th Mie-theory, but also a useful 
paramet r izat ion of the forward peak by reducing four parameters 
( character i z i ng , the size distribution) namely N, a»b, and y to 
only two, No and b. No is related to the total number of particles 
and b is related to the position of the maximum in number-density 
of this distribution. Hence, in our inversion scheme only two 
parameters appear, Nq and b, from which only the latter is needed 
for the knowledge of the forward peak, the most Important part of 
aerosol particle scattering. However, for longer wavelengths, 
when k is rather small, it is easy to find cases where the validity 
of the Born approximation applies throughout the entire range of 
0 ( 0 s 9 s 180°) and thus the formula (2-3) gives P](P) accurately 
for all scattering angles beyond the forward peak. 

It is easy to see from De i rmendj ian 1 s tables that this may 
not be true for shorter wavelengths. Then the method of separation 
(truncation by Potter (ref. 19 )) of the forward peak can be used. 

This is based on the favorable circumstance that the phase matrix 
P appears 1 i near 1 v in the equation of radiative transfer and that 
therefore it can be split into two parts (or more, if needed). 

Let Pf be the part of the phase matrix containing the "forward peak" 
elements, given by the analytical expressions (such as (2-3)) and P w 
the remaining part of the phase matrix such that P = Pf + P w . Since 
the higher order terms in the Lengendre series for the matrix elements 
needed for proper representation of the steep forward peak are in- 
cluded in the first part Pf, the elements of P w will then be repre- 
sented by a' Legendre series of fewer terms, or they can even be 
neglected with respect to the magnitude of the forward peak (see 
Fig. (1-2). 

In addition, we can exploit another advantage based on the 
linearity of equations of radiative transfer. They can be split 
into two equations, one containing terms related to the forward 
peak's part of the matrix, Pf, and the other to the remaining terms. 
Physically, the first equation defines the intensity and polariza- 
tion parameters as observed after scattering in the forward peak 
only . This is the well-known aureole that an observer sees as a 
bright ring around the sun. The equation for the aureole, exten- 
sively studied and described by Deirmendjian ( (ref. 11) p. lOOff), 
can be easily solved with the expression (2-3). The second equation 
for the radiation field in directions off the forward direction 
that remains after separating the equation for the aureole, retains 
the form of the original equation, except for a new term that 
represents the irradiation by the aureole. This suggests that in 
radiative transfer in a turbid atmosphere two sources appear: The 

sun ( 1 0° times brighter than the diffuse skylight) and the aurole 
(in some cases up to 1CP times birghter than the skylight). 
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Although this equation will contain the phase matrix with elements 
which can be developed in much fewer terms, it is still mathematic- 
ally very difficult to solve. Hence, the question arises of wheather 
it is necessary to solve these equations with different parameters 
for a "model library," or simply to replace the side-angle scattering 
by Rayleigh scattering with the same volume scattering coefficient. 

In such case only slight modification is needed of the computational 
scheme of molecular scattering based on the already existing tables. 

To compare these two methods, and to assess the accuracy of such 
an approximation is a considerable task. The time available in 
this contract has allowed us only to outline the basic principles 
of such an approximation. 

2.3 Effect of Ground Reflection . - In order to avoid unnec- 
essary complication, the assumption has been made in the previous 
text that the sky radiation transmitted by the atmosphere is com- 
pletely absorbed by the earth's surface and hence is not reflected 
upwards. Unfortunately, this happens only under rather rare con- 
ditions and the effect of so-called ground reflections has to be 
taken into consideration and proper modification established. 

For this purpose we need first to know the reflecting proper- 
ties of the ground. These have been measured from an airborne 
platform, primarily the total amount, of reflected energy from all 
directions - albedo or reflectivity. If the ground reflects accord- 
ing to Lambert's law, i.e.» if the reflected radiation is isotropic 
(equal in all directions) and unpolarized,, irrespective of the 
variation in intensity and the polarization of the incident radiation, 
such an albedo is all that is needed. For a more general case, it 
is necesary to introduce so-called bi-directional reflectivity, 
R(TT,Tr'), that gives the reflected intensity in the direction TT in 
relation to the intensity of the incident radiation from "FT'. If we 
consider also the state of polarization of the incident and reflected 
radiation, then R(tt, IT') has to be replaced by a 4 x 4 reflection 
matrix, in order to relate the polarization parameters of these two 
radiations. All these quantities have been measured by several 
investigators (ref. 20 and 21). The laboratory measurements indicated 
surprisingly similarity between different soil samples, and suggest 
that the ground reflection may be classified in several types that 
would considerably decrease its variability and facilitate simple 
parametr izat ion. In radiative transfer theory the ground reflection 
can be introduced rather simply, as it was first done by Chandrasekhar 
((ref. 1) p. 270 ff) as an additional source irradiating the atmos- 
phere from below. Experimental studies of bi-directional reflectance 
were first started by Coulson (ref. 22 ), extended for polarization 
by Coulson (ref. 23 ) snd Chen and Rao (ref. 21 ). Systematic studies 
of ground reflection in laboratories, and from airborne platforms 
have been conducted for several years by the UCLA radiation group 
under Professor Sekera's and Dr. Rao's direction with grants from 
the National Science Foundation, and a comprehensive article is 
being prepared for Advances i n Opt i cs . 
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Before closing this part, it should be mentioned that although 
the ground reflection adds a few new parameters to the "model library," 
these parameters differ in one point substantially from those related 
to atmospheric turbidity; except for seasonal changes, they are 
constant , and a few measurements 'under extremely clear conditions 
will be needed for their determination. 

2 . k Evidence for the Validity of the Proposed Inversion Method. 

A cursory reader may question the possibility of- determi n i ng from 
satellite observation the optical characteristics of the atmosphere. 

One must be aware that these observations are a consequence of 
complicated interaction of matter and radiation. In reality, satel- 
lite observations deal mainly with the reflection properties of the 
atmosphere. Ground based observations provide complementary infor- 
mation such as transmission properties of the atmosphere. 

Chandrasekhar's theory of radiative transfer (ref. 1) and subse- 
quent computations by other workers (ref. 20 ) established quantitative 
predictions of the diffusely reflected and transmitted light in a 
Rayleigh atmosphere. One reason for the deviations of experimental 
measurements from the 'theoret i cal predictions should be attributed 
to scattering by the particulates in a turbid atmosphere. Quite 
extensive skylight polarization measurements by the UCLA Radiation 
Group and others provided definite proofs that this was exactly so 
and that the observed deviations from the theoretical polarization 
of a molecular, Rayleight atmosphere are the most sensitive indica- 
tors of atmospheric turbidity. These classical, almost entirely 
qualitative determinations gave rise to the idea of changing them 
into quantitative ones, allowing estimates of turbidity parameters 
first from ground based and finally airborne or satellite observing 
platforms. Hence one should first look for the verification of the 
theoretical "inversion" method in ground observations. There are 
few such efforts, but one deserves mention in connection with the 
"forward peak" separation discussed above. Shortly after the 
polarization computations for a pure molecular atmosphere, and a 
dependable photopol ar imeter became available, Dr. E. De Bary (ref. 

25) performed an extensive series of skylight polarization measure- 
ments at UCLA with the intention of comparing the measurements of 
the polarization in the sun's vertical with the theoretical values. 

She found in a large number of days with moderate turbidity the 
curves of measured degree of polarization ran close to the theoret- 
ical curves, and substantial deviation appeared only near the sun. 

This fully supported the previously mentioned idea that the aerosol 
scattering in a turbid atmosphere can be well approx ina ted by 
forward peak only, considering the off side scattering as Rayleigh 
scattering with properly modified scattering coefficient. 

As another attempt, even more closely related to verification 
of the mentioned theoretical approach one should mention the work 
of Kevin Pang (ref. 25) who has applied the suggested method (ref. 

10) of estimating the optical thickness of Rayleigh atmosphere from 
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satellite measurement to skylight measurements. Results of his 
measurements proved completely the correctness of this method and 
the theories on which it was based. 

Finally, we should stress that careful measurements of the 
sun's aureole from platforms at different elevations are needed, 
not only to prove our forward peak analysis, but also to study the 
physically plausible invariance of particle size distribution in 
normal turbid atmosphere. Although it may seem that such analysis 
is irrelevant to the immediate goals of the contract, its value 
as a supporting study cannot be underestimated. 

Under another program, sponsored by the National Science 
Foundation, several attempts were made to measure the polarization 
of underlying atmosphere and the disappearing effect of surface 
features from high altitude balloons. These polarization measure- 
ments made from elevation of 92,000 to 95>000 feet over the White 
Sand Desert resemble much more those from a satellite since the 
thickness and hence the effect of the atmosphere above the obser- 
vation level can be neglected. Drs. C. R. N. Rao and T. Takashima 
undertook the difficult task of evaluating the amount of the aerosol 
particles and the Lambert reflectivity for several pairs of these 
quantities by matching the computation and observation in the sense 
of a "model library." This has been done under two simplifying 
assumptions, namely the "homogeneity" of turbidity (i.e., invariance 
of the ratio of the number of aerosol particles and the air molecules) 
and of the Lambert reflection. Both assumptions are to a great 
extent justified over a heated uniform desert. In Fig. (2-5) their 
results are presented in four wavelengths, the curves representing 
the computed distribution of degree of polarization in the sun's 
vertical, the dots are measured values. On each diagram two values 
are given, the so-called "turbidity factor" T=(?a + Tr VtL from 
which the total amount of particles can be derived, and the reflec- 
tivity A of the ground. The agreement of the computed and observed 
data are excellent not only in individual points, but also with 
respect to the identical run of the curves. A more convincing 
proof can hardly be found than that by a "model library" it is 
possible to find the fundamental parameters of a turbid atmosphere 
from satellite polarization measurements. 

The proof of feasibility to determine turbidity parameters from 
polarization measurements performed at lower levels brings with-it 
one difficulty, namely the knowledge of the Rayleigh radiation field 
within the atmosphere. The project director has studied the method 
of computation of the internal field of Rayleigh atmosphere and t 
the corresponding paper, entitled "Evaluation of the Internal Field 
in a Plane Parallel Atmosphere from Diffuse Reflection and Trans- 
mission" submitted for publication to the Journal of Quantitative 
Spectroscopy and Radiative Transfer. T. Takashima worked out the 
computational programs and computed the degree of polarization for 
the elements of the Convair 990 experiment described in Chapter 3-t 
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Computed values in Fig. (3-15) were plotted in dashed curves. Here 
again the reader must admit very good agreement with the measurements 
on a day with negligible atmospheric turbidity. 

2. 5 Other Applications of Polarization Measurements from a 
Satel 1 ? te . - In addition to many applications where the same quan- 
titative measure of atmospheric turbidity is needed (for example, 
transparency of the atmosphere for earth resources studies, effect 
of’ aerosol particles on temperature and humidity measurement (ref. 

26), determination of ozone, etc.) there are problems when the 
knowledge of polarization may be successfully used. In all these 
applications it is not. the diffuse reflection that is used for the 
analysis, but rather direct reflection with respect to the attenu- 
ation in the atmosphere. 

As an example one can mention the reflection of sun's rays 
over large bodies of water (glitter) where the inclination of the 
reflecting wavelets can be easily determined from Fresnel's law, 
if the degree of polarization of the reflected ray is known. Since 
in this law the refractive index of the liquid enters in a signifi- 
cant way, the polarization studies can be extremely useful in de- 
tection of oil slicks and simi 1 ar. surface features. All -these 
problems were discussed by several authors and will be summarized 
in the already mentioned comprehensive article being prepared for 
Advances i n Opt i cs . 

The law of reflection on the crystalline features of ice 
particles is more complicated. This will produce polarization of 
moderate values allowing a very simple distinction between layers 
of ice and water clouds. Theoretical and experimental studies 
leading toward such distinction are in initial stages at the UCLA 
Radiation Laboratory as part of a grant from the National Science 
Foundation. 

The application of polarization measurements from a satellite, 
mentioned above, are mere examples and do not exhaust all possible 
cases. Measurements of the polarization of light (of all wavelengths, 
from UV to IR) together with its intensity, will without any doubt 
be a necessity of the future, and judging from the overall and quite 
extensive interest in pol ar izat ion, probl errs by a series of Russian 
scientists dat i ng from- the last century (as may be seen from 
Rozenberg's (ref. 27 ) bibliography) we shouldn't be surprised if 
the new types of Russian meteorological or geophysical satellites 
were equipped with, similar photopo I ar imeters as described in 
Chapter 3* 
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3. THE EXPERIMENTAL PROGRAM 


Description of Equipment and Experimental Flights 


3.1.1 The Photopol ar imeter and Associated Equipment . - To 
facilitate discussion of the practical aspect of gathering data 
and calibrating the instrument we present here a brief operational 
description of the system. 


Operat i ona 
are the so-ca 1 1 
are defined in 
through an anal 
angle c o from t 
plate, introduc 
horizontal osci 
of light emerg i 
press i on : 

K cp)= \ [i 


1 Principles - The quantities that are to be measured 
ed STOKES polarization parameters. These parameters 
the following way: If a polarized light beam passes 

yzer with the transmission plane deviated by the 
he vertical direction and then through a retardation 
ing a phase difference 9 between the vertical and 
llations of the electric vector, then the intensity 
ng from the retardation plate is given by the ex- 


+ Q cos2ep + (U cos0 - V sin 0) sin2cp 


where I, Q, U, V are the STOKES parameters of the measured light 
stream. They completely define the state of polarization of the 
light stream, I = I, + l r being the total intensity, Q = \ f - l r > where 
1 1 , | are the i ntehs i ty components in the vertical and horizontal 

dfrect r ion, respectively. The inclination Y of the plane of 
polarization (containing the maximum intensity component) from the 
vertical direction is given by Q and U, ( tan2Y = U/Q). Finally V 
defines the ellipticity (V/l = sin2fi , where P = b/a, the ratio of 
the minor and major semi-axis of the ellipse described by the end 
part of the electric vector). If the ellipticity is small and 
can be disregarded, as is the case in atmospheric scattering, 
the retardation plate can be omitted, and the state of linear 
polarization is given only by three parameters, I, Q, U. Then as 
it follows from Eq. (3-1) for 0 = 0, the measurements of I for 

three positions of the analyzer are needed for the determination 
of these parameters. This can be simplified even more if the 
measurement in one position is replaced by the direct measurement 
of the intensity (without analyzer). Such procedure has been 
applied in a new photoelectric polarimeter (PPM) proposed to be 
used in this experiment. The basic parts are three barrels (see 
Fig. (3-1)) one (barrel i containing only a diaphragm filter of 
50% transmission, the other two containing an analyzer in a hori- 
zontal direction (barrel B) and an analyzer inclined by 45° (barrel 
D). Substituting the corresponding value ofm (90° for B and -45° 
for D) in Eq. (3-1) we see that the parameters Q and U are obtained 
by subtracting the intensities from barrels B and D from the inten- 
sity in barrel A. The subtraction of these intensities is being 
done electronically. For this purpose, a barrel selector, placed 
behind the analyzer, has an opening arranged in such a way that, 
as it rotates, it lets consecutively pass through the light from 
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barrels A, B, A, D, A, respectively. By a proper adjustment, it 
can then be arranged that the voltages are measured directly pro- 
portional to A-B and A-D, i.e., to the STOKES parameters Q, and U. 

The light beams from the three barrels are collected, after the 
barrel selector, by fiber optics light guides that illuminate the 
same photodetector after a passage through narrow band filters in 
a filter wheel rotated by a Geneva mechanism. The speed of the 
rotation of the barrel selector and filter wheel is' adjusted to 
give the measurements for all four filters in one second. In the 
fabricated model by Space Science Laboratory of General Electric Co., 
the fiber optics bundle has been bent back and the photodetector 
assembly mounted parallel to the barrels in order to reduce the 
length of the instrument. 

Subsequent to delivery on July 17* 1971* it was discovered 
that the polarizing prisms were installed 90° from their prescribed 
positions so that ro is 0° for B and -45° for D. This may be 
accounted for simply by making the inversion Q to -Q and U to -U. 
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Table 3-1 


Redundant Notation 


Horizontal equivalents in the table below are used 
interchangeably throughout this report. 


Old 

New 

Quant i ty 

I 

S 0 

Intens i ty 

Q 

S 1 



I 

Stokes Parameters 

u 

S 2 


q 

b 0 

Output voltages from 

U 

b l 

photopo lar imeter 

i 

b 2 


Q/l 

Til 


U/I 


Coordinates in 
T]-plane 

q/I 



u/I 



f 

f l 

Cal ibrat ion (trim) 
factors 

9 

f 2 
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Data Acquisition and Recording System . - The information 
from the photopol ar imeter goes to the mul t i p 1 . exi ng system, which 
arranges it serially together with a voltage indicating the scan 
angle and a voltage indicating the temperature of the optical mod- 
ule. The signals are then converted to digital form and recorded 
on tape. 

The beginning of a data cycle is indicated by a pair of pulses 
(referred to as the color pulses), one of which is generated by the 
position of the filter wheel and its "echo" generated from the 
first electronically. After the color pulses, the sequence is as 
f o 1 1 ows : 

qj, gain, scan, temp, 

Uj, gain, scan, temp, 
i^, gain, scan, temp, 
q 2 » gain, scan, temp, 

U 2 > gain, etc. 

The color channels are numbered from one to four in order of 
increasing wavelength, q, u, and i denote the values obtained 
directly from the instrument. The actual Stokes parameters (Q, U, 
i) are calculated by taking into account the measured "zero" levels 
and the instrumental trim factors (these latter quantities are 
determined by independent calibration). 

Thus each complete cycle (1 second) contains 48 pieces of 
information. After 15 cycles, the time is recorded and a new block 
begins (The total block period is thus 16 seconds.) The basic 
computer program averages the scan angle and temperature for each 
Q, U, I sequence and calculates the corresponding degree of polar- 
ization and plane of polarization. 

Calibration Factors - Formalisms . - In order to discuss the 
quantities that characterized the actual behavior of the photo- 
polarimeter, it is necessary to introduce an extended notation. 
Results of attempts to measure the quantities discussed will be 
given in Sect. 3.1 . 2 . The color channel will be supressed in this 
section, the results being applicable to all colors. 

The Stokes parameters for the incident radiation are desig- 
nated by Sq, s,, s 2 > and s,, corresponding to I, Q, U, and V in 
the previous notation. The output voltages from the three barrels 
will be designated b n , b,, and b 0 , corresponding to I, 2B, and 2D 
(see Table 3-1 ). 0 1 2 

In this program it will always be assumed that s, = V = 0, so 
we will deal with the three component vectors * 
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where U is a 3 x 3 matrix defined as the instrumental matrix. 
Ideally, U would have the form 



(3) 


according to the instrument design principles. (We have taken 
into account the 90° mi sor i entat i on of the prisms here. (Note: 

In all that follows, the calibration of the I barrel is taken as 
unity. ) 

Since the three barrels do not have the same transmission 
characteristics, the ratio of transmission of the b, or b ? barrels 
to the b« barrel must be taken into account. These are designated 
fj and f^> respectively, and U then has the form 



The processor performs an operation represented by a matrix tt 
where the output v if given by 


v 


=nb = ttUs 


(5a) 
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where i deal 1 y. 


r 


0 0 
1 0 
0 1 


(5b) 


The calibration problem is to find the calibration matrix C 
such that 


C = (mU)" 1 = K" 1 


( 6 ) 


If rr and U have the ideal forms given by Eqs. (3-4) and, (3-5) we 
have only two factors to determine for each color channel. K will 
have the form 



it is useful to define tt by 


_ s 

1 = ¥ 0 (9) 

Since TU = 1 always, it is convenient to consider as a two- 

component vector. The m- plane is the plane with and ^2 
as coordinates', T| being taken here as the abscissa! Recalling 
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that 



s 

0 


tan2Y = - — 
s 2 


( 10 ) 


where P is the degree of polarization (which we will call simply 
"polarization") and y is the angle of polarization (measured 
clockwise looking at the light source), we then have 


P 



2 

+ ^2 


:! i 

ta n 2y = ... (11) 

Tl 2 


Thus, in then -plane the vector from the origin to a given point 
has a magnitude equal to the polarization, and the angle with 
respect to the abscissa is twice the angle of polarization. Thus 
P and 2Y are polar coordinates corresponding to the Cartesian 
coordinates ni and T>2 (Figs* (3-2) & (-3). We now have three 
equivalent forms for expressing the state of polarization of the 
1 i ght : 


r 



V s , * 

s 2 

_ . 



1 1 ^ 


<P 


2 y 

v V 


For our purposes here, the latter two forms describing the 
behav ior of the light in theTl -plane will be most useful for d is 
cussing the calibration and behavior of the instrument. We have 
for example . - 
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cos2f 


T| = P 


s in2y 


(= P) 


( 12 ) 


Consider the vector v defined analogously to Ti by 
v 0 

and using Eqs. (3-5a)and (3-8) 
v = K + M 


(13) 


(14a) 


where 



(14b) 


(14c) 


It is convenient to maintain the same symbols for certain 
matrix operations in reducing to the Ti-plane description. There 
should be no confusion since d i f ferent symbols are used for the 
state vector in the T-plane. The requirements for the reduction 
from the 3-d Stokes space to the 2-d Ti-plane of an arbitrary matrix 


M = M. , (i,j = 0, 1, 2) 
1 » J 


is that 


M 0, 1 = M 0, 2 - 0 


(15a) 
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Then the 3X3 matrix operation is replaced by a 2 X 2 matrix 
operat i on. 


M • M. , 

i, J 


( i,j . 1,2) 

'M, 


plus a translation m =1 


i . e. , s' = Ms 


M 


20 / 


(15b) 


( 15c) 


becomes 


V = Mil + m. 


(16) 


For successive operations, 


s' = NMs 


( 17) 


on the Stokes vector we have 


T»* = T| + N(m + Mm) 


= 71 + Nm + NMT) 


(18) 


on the 71-vectors. As an example, the product, matr ix CK from (3-7) 
and 3^8) becomes . 

/’ °\ 

CK — » c + C(k + K) = c- c + CK = 

\ 0 1 


/ 


09) 


(Note that if we had designated C as K ^ we might be tempted to 
reduce 


K " 1 s« 


-T< + K " T| 


( 20 ) 
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which would be incorrect.) (The validity of the assumption 
eq. (3-1 5a ) will be discussed later. This implies that the 
reading from the l-barrel is invariant to the polarization of 
the incoming light.) 


Returning to Eq. (3~1*0> the form of K given is not general. 
We choose to write the general form as 



f l T 12 


f 2 T 22> 


(21 ) 


which is consistent with our definition of the f's. We will call 


Tf = the state vector or "state" (or polarization vector); 
v = k + Klf = the output vector or "output"; 

H* = c' + C'v = the computed state vector or "computed state", 


in order to distinguish the various vectors we will be dealing 
with. The last quantity defines a situation in which we are using 
an incorrect calibration matrix constructed analogous to C but 
using incorrect calibration factors. 

Special Cases . - If the instrument is rotated about its axis 
while viewing a source whose polarisation remains fixed during the 
time of rotation^, the state vector T) will move along a circle of 
radius P in the T|-pl'ane centered at the origin. For a 180° rota- 
tion of the instrument, T) will make one full transit of. the 
circle and for a full 360° rotation of the instrument T| will travel 
twice about the circle. 

If the instrumental matrix is not unity, the output vector 
will describe some different closed curve. (In full generalities 
the curve is. required to close on itself only after two full ro- 
tations of T|. ) We will now discuss how certain special cases affect 
what happens to the circle. 

Case 1 : K 12. = K 21 7 °‘ The expected case with no coordinate 

mixing but unequal barrel transmissions. From eqs. ( 3 ~ 1 *0 and 
(3-12) we see that the curve generated by the output vector 
will have its center displaced by the vector K. There will also 
be an elliptical distortion of the circle. 

The equation for this ellipse is 
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2 


( 22 ) 




Note that it is not necessary to determine the elliptical properties 
of the curve in order to determine f] and f 2 » since these quantities 
are fully determined by locating the center. This suggests two 
methods of determining the f's: 


a) 90° Intervals - Taking a measurement at an arbitrary 
setting, rotating .the instrument through 90° and taking another 
measurement. Then 


b ( Tq ) + b(v Q + 90°) = 
2 



(23) 


b) Rinq-A-Round - (RAR) - Rotating the instrument through 
360° continuously and graphically determining the center of gravity 
of the resultant data curve. 

While these two methods are ostensibly the same, the first 
requires an accurate knowledge of the angle of rotation of the 
instrument. A variation of the first method was investigated 
in order to attempt to make use of in-flight data in determining 
the f's. This method which follows requires knowledge of the di- 
rection of polarization (y ) but not the degree (P). 

c) , Spec i a 1 angl es - When y =+90°, 0° then from (3-12) and 
( 3- 1 4) we f i nd 

iMO 0 )) 

\=f 2 -l 

b 2 (90°j 

and if Y =+45°, (24a) 


^(+45°) = f,-l. 


(24b) ' 
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This method requires an accurate knowledge of Y since an error 
6Y will contribute an error P5f to. the measurement ( si n SY 6YA 
On the other hand, since cos6Y «1 - illi— , a more reliable value 
might be obtained from 2 

b,(0°) + bj (90°) . fi V (25«) 

2 


and 


b 2 (45°) + b 2 (-45°) _ f„-i. 

2 z 


(25b) 


If. the angles are slightly off from the desired values but known, 
we would have to solve, two simultaneous linear equations: 


b i ( Y ' ) = (fj-1) + f|P cos2y' 
b 1 ( Y ' ' ) = ( f 1 - 1 ) + fjP cos2Y 1 


(26a) 


for fj and P, and 


b 2 ( Y 1 ) = (f 2 -l) + f 2 P sin2y' 
b 2 ( Y" ) = (f 2 -l ) + f 2 P sin2y' 1 


(26b) 


and f 2 and P. 

As before, we would expect a better value for f, ify 1 
and Y 1 1 are near 0° or 90° and a better val ue for f 2 i f Y 1 and 
Y 11 are near ±45°, since the effects of errors in' the 1 angles' would 
then be minimi zed. 


Case I I! - Effects of prims misalignment (mi sor i entat ion) . We 
will show by a series of transformations that the general result of 
prism misalignment may be considered as a rotation in the Tl -plane 
plus an elliptical deformation. 

The essential axes in real space are defined in Fig. (3-2). 

(T) | ) is primary reference plane of the instrument corresponding to 
Y = 0 plane of polarization. We label it (n 1 ) since the correspon- 
ding vector i n T| - space will be along the Tl ^ - Px is. The prism of 
the B-barrel should have its transmission plane or i ented a 1 ong 



(n,)» but if it is misaligned we will label its error angle by 6 
as shown. The same statements apply for the axis (T| „ ) and 6 is, 
the error angle of the D prism. L 

The axis (C , ) is the b i sector of (ti . ) and (t| 2 ) » and (q , ) is 
the bisector of (B) and (D). In the T] -p 1 ane, £ , Ts located at 2y 
= 45° and its conterpart ' s at .= 135°* 

If the f's have been properly accounted for, we may assume 
them to be unity in what follows. 


From Eq. (3-1) we see that the effect of error angles e and 5 
on the output voltages is given by 


b = 


0 

cos2e 
-s in2 6 


0 

s i n2 6 
cos2 6 


or in theTi -plane by (3-1 4a) with 

/ . \ 

/ cos l b s 1 n2 h \ 

-■'■1 

\-s i n2 6 


cos2 6 


and k = 0. 


(27) 


( 28 ) 


If we reduce the data without properly accounting for T, then 
we will be plotting 


Tl' = T'n, (29) 

i.e., our plot axes Tlj, t| 1 , will not be identical With T|j, T) , 
as desired and the curve ^generated in a rotary-type experiment 
will no longer be a circle. To see what effect T has, we first 
define / 


e 


is + e ; „ _ 6-6 

2 Y " 2 


(30a) 


and the inverse transformation, 


6 = 


e + y; 8 = e - y 


(30b) 




and substitute Eq. (3-30b) into Eq. (3-28). Define the rotation 
matrix A(0) by v 


A ( 0 ) = 


cos0 sin0 
l-sin© cos0 


we notice that T can be factored into 


T = G(2 Y )A(2 e ) 


where 

G(2 y ) 


•r 

cos2 y 

-s i n2 Y 


-s i n2 Y 
cos2 Y 


which we will call the deformation matrix. 


A( 2? ) is the rotation mentioned earlier; it rema 
the character of the deformation. We note that G(2 Y ) 
by the transformation / , 


G 1 ( 2 y) = A(J)G(2v)A(-{) = 


cos2 Y s in2y 

0 cos2y 




corresponding to the transformation 

C = A(J)tJ' = A({)G(2 Y )A(-g)A(J)A(2 e )Ti 


Noting the inverse 


G( 2 Y ) _1 


1 

Det G' 


G ' ( -2 y) 


cos2 Y + s i n2 Y 


cos2 Y 


(3D 


(32) 


(33) 


ns to show 
is diagonal ized 

\ 

0 

s i n2 y 

(34ab) 


(35) 


(36) 

0 

J 

- s'in2 Y 


0 


l 
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or 


G'Uy)" 1 C = A(?)A(2e)T|. 


(37) 


Taking the magnitude of the vectors on both sides we have 
£2 Q2 

I 2 2 

( cos2y + s i n2y)^ + (cos2 - sin2y)^ ^ J 


(38) 


recalling that the matrix A( 0) does not alter the length of a 
vector. We now see clearly that the circle is distorted into an 
ellipse with the ratio of major to minor axis given by 


a_ cos2y + s i n2y 
e - b - cos2 y - s i n2y 


\ . 

1 + 4y for ( y « 1 ) . 


(39) 


Rotating to the £' system, 

A ( 1+ ) iT = A(-2G)G'(2y)" 1 -^- 


(40) 


and finally 


n = 


A( -5)A(-2 e )G'(2y)“ 1 


A(7t ) T| 


or 


c = A(4)A(-2G)G'(2y)" 1 A(J) 


(41a) 


(41b) 


giving the calibration matrix in terms of matrices already ex 
pi i c i t 1 y stated. 
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Case III : Non-specific component mixing. The case of a 

completely general form for K is applicable to the instrument in 
the scan platform housing, where the gross effect of the quartz 
window including reflections , and the instrumental optics must be 
accounted for. In general we would expect 


K,2 « K n ; K 21 « K 22 . 


From the form of Eq. (3-1^+a) we see that even in the most general 
case, method (b) under Case J_ is appl icable for determini ng 1<. 
Accounting then for f] and f 2 > we have to 'determ'i ne the matrix 
T. While measurements at two known angles (y) could determine the 

ratios T 1 2 /T 1 1 » T 2 1 ^ T 2 2 ’ T 1 1 2 2 ' we wou1d have to know the de 9 ree 
of polarization in order to determine the individual values. In 
fact, if we take measurements at three known angles, the quantities 
which we can determine directly are 


-K 

v i ? i » 


f 1 PT 


IT f 1 PT 1 2 

•r nr 

21 ’ ' V ' 22 . 


The statement above follows from this. 

Effect of Errors ? n the Cal ibrat ion Mat r i ces . - Ifthecali- 
bration factors used for computation are incorrect we will have 

T} = c 1 + C' v = c' + C ' K + C' K 71 (42) 


The center of the closed curve described as the instrument rotates 
about its axis is displaced from the origin by the vector 


e = c' + C'k (43) 
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If the difference is sufficiently small, we can write this as 



\ 

, 4 


(44b) 


This leads to an uncertainty in P of ±lel . Assumi ng ! that K has 
the form of Eq. (3-7) letting f-f] = 6f» we have 



1 + 



(45) 


leading to an error in P of 


6P 

P~ 


1 2 

—1 cos 4 2y + 
T 1 



sin 2y. 


(46) 
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Note that this error is proportional to P, while the effect above 
was not. Since P is always less than 1, the e term will always 
dominate the error, and for low polarization ( < 0.2) we can ignore 
this distortion altogether. 

It is then possible to interpret the effect of this error 
by geometric interpretations in the T] -plane (see Fig. (3"3))* As 
an immediate example, if we measure a value of P to be 0.10, but 
do not know our f's to hott er than 1%, then we have an uncertainty 
of +14% in P (6 P = e ) or P = 0.10 +0.014. 

3.1.2 Aircraft Installation . - Purpose of the Polarimeter 
Scan Platform - The scan platform serves a twofold purpose: (T) 

By facilitating the scanning motion of the polarimeter, it makes 
possible the taking of quas i -s imul taneous measurements of the 
radiation field in several directions; (2) It provides means for 
aiming the polarimeter continuously on selected ground targets 
during the flight over the designated territory. Consequently, the 
important contribution of the ground reflection to the radiation 
field in the atmosphere remains an invariant function as seen by 
the polarimeter during the scan over the target area. 

Furthermore, the scan platform (with the other aircraft 
instrumentation) provides means for a check on the accuracy of 
the determination of the polarimeter I i ne-of-s i ght nadir and 
azimuth angles at different local sun elevations. 

Scan Platform Design and Location of the CV-990 . - The 
polarimeter scan platform described below was designed and con- 
structed as an interface between the UCLA skylight polarimeter 
with a bore-sighted television camera and the experimental aircraft 
Convair 990. The design requirements were postulated as follows: 

1) Full 360° scan capacity. 

2) Means for aiming the polarimeter on a target and TV re- 
cording of the target environment. 

3) Provision of a pressurized and thermally insulated 
compartment for the polarimeter and the bore-sighted TV 
camera. 

4) Streamlined shape of the exposed section to prevent any 
vibration due to aerodynamic forces. 

5) Ability to withstand the aerodynamic forces at extreme 
flight-load factors of the aircraft 

6) Electronic remote control of all functions of the system. 

7) Easy access to all parts for installation and calibration 
purposes. 
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In a previous experiment (ref. 28) (Grant No. 05-007-041) a 
polarimeter was located at the very end of the tail section of the 
aircraft. However, this arrangement did not satisfy requirement 
number 1, and because of considerable vibration of the aircraft 
structure in the tail section, requirement number 4 was not complied 
with either. 

Several design arrangements were studied, proposed, and dis- 
cussed with the NASA/Ames Research Staff. The design presented 
schematically in Fig. (3-4) was finally accepted, and the construc- 
tion executed partially at UCLA and at NASA/Ames Research Center. 

The scan platform is installed on a "low-boy" supporting structure 
in the front part of the passenger cabin. It consists of an 
aluminum (6061 -1'6) alloy tube, 6 inches in diameter and 5 feet long, 
mounted between two special ball bearings on the supporting struc- 
ture. To one end of the tube is attached an enclosure containing 
the primary electronics for processing of the electronic signals 
from the polarimeter. On the same end is also mounted a large 
gear connected by a link chain to a small sprocket of the scan 
drive motor (Fig. (3-5a))» An angular resolver (5 revolution pot.) 
is connected to the sprocket shaft. The scan drive motor is 
mounted on a spring-loaded plate, sliding between two guiding rails 
(see Fig. (3 - 5b)) in order to eliminate any slack of the link-chain. 
The other end of the six- inch tube penetrates through a side window 
of the passenger cabin and holds by means of a flange, a rotation- 
ally symmetrical streamlined fairing. A rotational seal, axially 
sliding on the six-inch tube, is pressed by the cabin air against 
the cabin window adapter plate. The polarimeter with bore-sighted 
TV camera is mounted inside a sealed and electrically heated com- 
partment centrally located in the fairing. Four optical windows 
(one for each barrel of the polarimeter and one for the TV camera) 
sealed in the wall of the compartment transmit the radiation to 
the instruments. The shape and relative size (30 inches in diameter) 
of the fairing are apparent from Fig. (3-6)-. Another small, stat- 
ionary fairing is mounted to the aircraft hull in close proximity 
to and downstream of the exposed part of the outside penetrating 
six-inch tube (Fig. (3-4)). As is apparent from Figs. (3-4) and 
(3-6), the scanning motion is in a vertical plane parallel to the 
center line of the hull of the aircraft. 

Although the scan platform was explicitly designed for use on 
the Convair 990, its lightness makes it applicable for use on a 
smaller aircraft, helicopters, and balloons as well. 

Scan Platform Control System and Instrumen t ation . - The 
scanning motion is provided by a motor -genera tor (CRS NC-100) and 
feedback servo-amplifier system (CRS NC-102). By means of a hand- 
operated potentiometer any scanning speed in forward or reverse 
direction can be set (constant and independent of the torque) be- 
tween 0.1 to 100 degrees per second. The monitoring instruments 
and the electronic data acquisition system are mounted on a "high- 
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boy" standard instrument rack (Fig. ( 3 ~7 ) ) • The scan position in- 
dicator has scale calibration marks for scan zenity and nadir 
angles and also for 60° backward from nadir positions. The, scan 
speed indicator has scale calibration marks for initial target scan 
speed, sky scan speed, and return speed values. Two limit switches 
automatically stop the scanning motion when in either direction the 
polarimeter approaches horizontal forward (flight-oriented) position 
(Fig. (3-8)). 

• i 

. The polarimeter 1 i ne-of-s i ght view is obtained by means of a 
bore-sighted TV camera (Sony ACV-3400) displayed on a TV monitor 
and simultaneously recorded on a v i deo- recorder (Sony AV-3400) for 
later reference (Fig. (3-7)). The full monitor screen shows the 
TV camera field of view, and the small circle in the center indicates 
the polarimeter field of view. The angle of view of the polarimeter 
is 3 > while that of the TV camera is 12°. Therefore, an area 
linearly four times larger than the target can be surveyed and 
possible interference by, for example, nearby drifting clouds, 
identified. For the "sky scan," a neutral density filter (Kodak 
ND-4) is placed in front of the TV camera lens (Fig. (3-9)). This 
permits scanning in the sun-vertical "through the sun" without 
damage to the Videcon tube. This feature also facilitates a 
double check for accurate determination of the heading of the 
aircraft. During the ground (target) scan, a rotary solenoid is 
energized to swing the neutral density filter out of the 1 i ght 
path. The sound track of the videotape serves for recording the 
comments by the flight coordinator or any member of the operating 
crew as well as for recording of time signals (WWV or from a separ- 
ate tape recorder) for later identification of the video record 
with the corresponding data output from the polarimeter. 

A rotary solenoid-operated shutter is mounted in front of 
the three optical barrels of the polarimeter for in-flight calibra- 
tion (zero signal) check and for the elimination of unusable data 
which would be acquired during "in between measurements" periods 
(Fig. (3-9)). Also visible in Fig. (3-9) is a miniature blower 
which circulates the air (mostly dry supplied continuously from 
storage bottles) inside the instrument compartment in order to 
maintain uniform temperature of the instruments and to prevent 
fogging of the optical windows. The instrument temperature is 
monitored by means of a thermistor (attached to the body of the TV 
camera) and recorded multiplexed with the data output from the 
polarimeter. The temperature level is controlled by four-step 
switching of the current to the heating stri ps attached to the 
walls inside the fairing compartment. 

A second auxiliary TV monitor is mounted on top of the instru- 
ment rack (Fig. (3-7)) to display the view of the target area as 
seen by a TV camera located in the nose of the aircraft. This 
camera is equipped with a remotely controlled zoom lens, and the 
displayed picture is helpful for identifying the target prior to 
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the initiation of the ground scan. 

Testing and Calibration of the Aircraft installation . - Upon 
the assembly and installation of the scan platform, controls and 
data acquisition system on the CV-990, a test flight was performed 
to determine whether the scan platform was air-worthy and whether 
the fairing affected the performance of the aircraft. Several 
accelerometer and pressure pick-ups were placed inside the instru- 
ment compartment and on the inside wall of the aircraft cabin in a 
downs t ream.-.d i rect ion from the scan platform. Analysis of record- 
ings indicated appreciable vibration of the cabin wall, caused 
probably by. the turbulent flow between the fairing and the aircraft 
fuselage. To eliminate this problem, the NASA engineers designed 
and additional stationary fairing (Fig. (3-4)) mounted directly on 
the window plates and extending from the six-inch tube to the 
second window downstream of the scan platform. A second flight 
test yielded results indicating an acceptable flow pattern in the 
wake of the fairing. The scan platform and the instruments were 
entirely vibration-free, as was also indicated by a completely 
steady picture. on the TV monitor screen. 

The scan angle calibration was performed by means of a highly 
sensitive precision water level attached to the polarimeter at 
four positions of the instrument, as aiming at zenity, nadir, and 
head or taii horizon. The calibration curve showed excellent 
linearity and accuracy of readout to 0.1° (Fig. (3-10)). The 
polarimeter electronics and data acquisition system as well as the 
television circuit were tested for interference from the aircraft 
power system. No effect on the polarimeter electronics and only 
a very slight effect on the TV monitor picture were noticed. 
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Fig. ( 3 -5 a ) -- Photograph of the Interior Part of 
the Scan Platform Installation from the Cabin End. 



Fig. (3 
of the 
End of 


1*0 Photograph of the Interior Part 
Scan Platform Installation from the Aft 
the A i rcraft. 



Fig. (3-6) -- Exterior 
Platform Instal lation 


Photograph of the Scan 
on the CV-990. 



Fig. ( 3 - 7 ) -- Photograph of the Instrument 
Rack and Control Panels for the CV-990 
Instal lat ion. 
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Fig. (3-8) -- Schematic Defining Scan Angle and Regions of Scan. 


Fig. (3-9) -- Photograph of the PPM-TV Camera Bores ighted 
Package. Note the shutter in front of the camera lens containing 
an ND-4 filter to protect the vidicon tube. 




E (Deg°) 








3.1*3 Supporting Data Systems - Direct Sampling of Aerosol 
( NCAR) . - Direct sampling of the aerosol was carried out under the 
direction of Mr. I. Blifford of the National Center for Atmospheric 
Research), Boulder, Colorado. Since previous commitments prevented 
Mr. Blifford. from being on the flights, his equipment was operated 
by Mr. G. Dolan, also of NCAR.. 

V : 

The impactor used by Mr. Blifford has been described in the 
literature and its installation on the CV-990 presented no diffi- 
culties. (ref. 29 ) I n addition, the NCAR workers installed a 
filter collection system and an optical counter. The optical 
counter malfunctioned on the Carson Sink flight and despite Mr. 
Dolan's efforts, remained inoperable throughout the rest of the 
flights. This instrument would have given a real time indication 
of particle densities had it functioned properly. 

After the flight the filters and impactor slides were returned 
to NCAR for analysis. Th.e results are presented in Section 3*2.1. 

Aircraft Support Systems . - NASA/Ames provided numerous services 
and systems which were useful and/or indispensable in accomplishing 
the experimental program. This section does not present a complete 
list, but only those which were essential for data reduction and 
interpretation both during the CV-990 flights and afterward. 

1) ADDAS - The Airborne Digital Data Acquisition System on 
the CV-990 consists of a Hewlett Packard 2 1 1 6/ B/2 1 5 OB combination 
of a general purpose computer with extender plus periphery equip- 
ment. The system was utilized to record the output from the preamp 
as well as the scan voltage, etc., as a back-up data recording 
system for the experiment. The system records on a 9-track digital 
magnetic tape 800 bpi at 37*5 ips. Also recorded were the signifi- 
cant aircraft parameter outputs from the INS (described below). A 
real line printout of selected aircraft data was made on a Nortec 
200 lineprinter every 10 seconds (or during the r i ng-a- rounds , every 
second). The printouts were presented to us after each flight. 
Duplicates of the aircraft tapes covering the time period of data 
taking were also provided. 


2) INS - The Litton 51 Inertial Navigation System provided an 
update every 1.2 seconds of the following aircraft and navigational 
parameters: 


a) distance to go 

b) time to go 

c) cross back distance 

d) desired track 

e) track angle error 

f) drift angle 

g) align status 

h) present position latitude 

i) present position longitude 


J) 

ground speed 

k) 

track angle 

1) 

true 

head i ng 

m) 

wi nd 

speed 

n) 

wind 

durat i on 

o) 

long 

i tude waypo i nts 

P) 

1 at i 

tude waypoi nts 

q) 

rol 1 

angle 

r) 

pi tch angl e 
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3) Others - Ground-based computers at NASA/Ames w ere used 
to spot check the UCLA data tapes immediately after each flight. 

A belly camera recorded the terrain below during the flight. 

An on-board V i s i corder was used for periodic real-time checks on 
the PPM preamplifier to verify that the equipment was operating 
properly.' 

Other Experimenters . - In addition to the NCAR direct sampling 
equipment* other experiments were on board and run on some or all 
of the f 1 i ghts : 

1) Stratospheric Sampling Program - J. G. Poppoff, NASA-ARC. 

2) Laser True Air Speed Program - Ri M. Munoz, NASA-ARC. 

3) Infra-Red Water Vapor Measurements - P. M.' Kuhn, NOAA-APCL. 



3 • 1 • 4 The Convair 990 Flights - Review of the Flight Program . - 
The UCLA, team participated in a total of five flights on the Convair 
990 at NASA/Ames, fourcf which were prime flights for UCLA. The 
original schedule called for the flights to begin on November 16, 
but failure of one of the aircraft engines forced a delay in the 
schedule. Simultaneously, a malfunction of the tape recorder was 
detected and traced to an open winding in the head. The tape 
recorder was taken to the Kennedy Company for repair, then rein- 
stalled on the CV-990. A summary of the flights is as follows: 

December 13, Carson Sink, Nevada . - This site was chosen for 
a test site due to poor weather conditions and restrictions at 
preferred target areas. The first target site was abandoned when 
low-lying cumulus clouds drifted over the site. Scattered snow and 
high winds at the higher altitudes made it difficult to track a 
homogeneous target. 

After the flight it was found that tape errors had developed 
after about an hour of operation. The tape was taken to the Kennedy 
Company in Altadena for analysis and the following day the tape 
recorder and A/D converter were taken to Kennedy, but no instru- 
mental malfunction could be found. 

December 17, Canada. -'The equipment was reinstalled on the 
CV-990 and the UCLA experiment rode piggyback with a stratospheric 
sampling experiment in order to test equipment. A test signal 
box borrowed from Kennedy was taken along, and alternate recordings 
from the test box and data system were made. From this flight it 
was decided that the errors were resulting from overheating of the 
A.D converter circuits. Simultaneously, the programmer came up 
with a program to read the bad tapes, but running the program is 
not only expensive, but consumes about an hour of supervision for 
every minute of CPU time. 

December 21, Great Salt Lake Desert . - Again low level clouds 
moved in over the target area shortly after our arrival, but the 
patterns were continued to test equipment and routines. The A/D 
converter was turned off when not in use and its ventilation improved. 
A test also was made wherein the output of the photopol ar imeter 
processor unit was recorded on the aircraft analog recorder. On 
return to Moffett Field, the analog tape was played back through 
the remainder of the UCLA system and transformed to the usual digi- 
tal record. This method was successful, pampering the A/D converter 
was also successful and the other data tapes were error free. The 
ventilation was further improved by running a hose from the air- 
craft's air conditioning system directly to the A/D converter. 

December 29, Mexic o. - The weather conditions cleared unexpec- 
tedly and permitted a flight to the Mexican desert (Gran Desierto) 
just south of Yuma, Arizona. Again, high winds were encountered at 
40,000 feet but the target area was homogeneous over such a large 
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region that the resultant drift is not expected to be reflected in 
the data (see Fig. (3-11)). The weather was exceptionally clear 
with neither cirrus above nor any clouds on the horizon. The hori- 
zon view indicated light haziness and the presence of an inversion 
northward. ■ 

December 3 Q> Mexico. - Since weather conditions were much the 
same, the preceeding day's flight was repeated. Several high cir- 
rus clouds were present over the target area and there was an 
obvious increase in haziness though it was still light. No cirrus 
clouds were observed to obscure the sun. The terrain of the Gran 
Desierto appeared exceptionally flat and homogeneous for an area 
of about 400 square miles and was superior in this respect to any 
areas overflown in the United States. For this reason as well as 
because of the difficulties already related, only data from the 
Gran Desierto flights have been reduced. These two flights are 
referred to as Yuma I and Yuma II, respectively, or Y1 and Y2, in 
the brief notation listed in Table 3“2. 
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Table 3-2 


-- Key to Code or Brief Notation 


Symbol 
Yl, Y2 


Rl, R2 


LI, L2 , 
L3 , L4 

Cl, C2 
C3 , C4 


Des i gnat ? on 

flights over the Gran Desierto on 
12/29/71 and 12/30/71, respectively 
(Yuma I and Yuma II) 

rosettes at 10,000 ft and 40,000 ft 
altitudes, respectively 

legs 1, 2, 3, 4 at nominal 0°, 90°, 

135°, and -135° to the sun, respectively 

color channels centered nominally at 
3800 A, 4400 ft, 5000 A, and 5800 ft, 
respectively (See Table 2 in the Appendix). 
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Flight Experiment Procedures . - The flights prior to the Yuma 
flights were invaluable for debugging equipment and establishing 
suitable flight procedures. For example, it was found impossible 
for the scan operator to recognize and pinpoint the target area 
on the TV screen as the plane approached from different directions 
at different altitudes. Instead, for the Yuma flights the INS was 
utilized for target acquisition. 

The flight procedure sequence used was as follows: 

1) The plane flew over the area at 1,000 feet to select a 
suitable target. When the plane was directly over such a point, a 
button was pressed, locking this point into the INS as ground zero. 
During this time aerosol samples were collected and after "ground 
zero" selection, the plane flew a "racetrack" pattern in the plane 
of the sun until sampling was completed. 

2) The plane ascended to 5,000 feet for more aerosol sampling, 
again in a racetrack pattern. During the 1,000 and 5,000 feet 
patterns, the photopol ar imeter scanned the sky. 

3) The plane ascended to 10,000 feet and began the rosette 
pattern shown in Fig. (3-12). For legs 1, 2, 3, and 4 the desired 
relative solar bearings are 0°, 90°, 1 3S°» and -135°, respectively. 
On the approach to each leg, the scan operator set the instrument 
scan angle 30° below the horizon (60° forward of nadir). When the 
INS indicated the plane was the desired distance from ground zero, 
i.e.» when the instrument and TV were centered on ground zero, the 
Operations Officer relayed the signal to the scan operator who 
began scanning to keep fixed in the center of the TV screen what- 
ever was there at the time. When the INS indicated that the target 
was 150° below the horizon ( 60° back of nadir) the signal was again 
relayed and the scan was set to a fixed 2.5°/second rate, in order 
to scan the downward skylight radiation. Simultaneously the ND4 
filter was closed on the TV camera. The airplane continued a 
straight course until the sky scan was completed and then began the 
approach to the next leg. Between data legs the instrument shutter 
was closed to record zero levels. "Zerd 1 and "data" portions of 

the tape were separated by file gaps. (See Fig. (3-8)) 

b) After all four legs were completed, the same procedure 
was repeated at 40,000 feet, except that area restrictions forced 
curtailment of the sky scan on leg 4 at 270° (zenith) on 12/29 and 
at the horizon on 12/30. 

The Rinq-A-Rounds (RAR). - The only feasible method for deter- 
mining the calibration factors for the instrument as installed was 
to rotate the scan platform until the instrument was pointed at the 
zenith and have the CV-990 taxi in a circle on the hangar apron 
at NASA/Ames. 
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Ideally, a clear blue sky was desirable for this purpose but 
such conditions were never present over Moffett Field during the 
experiment period. The conditions for the r i ng-a- rounds (RAR) that 
were performed were: 

1) Salt Lake RAR - prior to flight, with uniform overcast 
overhead. 

2) Yuma I RAR - after flight, clear skies directly above but 
many scattered clouds immediately surrounding. 

3) Yuma II RAR - after flight, conditions similar to those 
for Yuma I RAR. 

The windows were at their cleanest for the Salt Lake RAR, 
having just been installed the night before. Prior to the Yuma I 
flight an attempt was made to • remove the more obvious dirt particles 
with a long cotton swab. No attempt at cleaning was made between 
the Yuma I and Yuma II flights. 

After the Yuma II flight, the ports of the scan housing were 
sealed to preserve the dirt to permit further calibration studies 
at UCLA. 
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Fig. (3-12) -- Rosette Patterns Executed at 10,000 ft. and 
40,000 ft. The Legs of the Rosette Are Defined with Respect 
to the Sun at 0°, 90°, 135° (-45°) and 45° for L 1 , L2, 

L3 , and L4, respectively. 



3. 2 Results of the Flight Experiments 


3.2.1 Direct Aerosol Sampling . - In order to have some contin- 
uity in the discussion of the reduction of the polarimeter data we 
present first the results of the direct aerosol sampling. At the 
conclusion of the flight programs the slides and filters from the 
d i rect -.sampl i ng experiment were analyzed in Mr. Blifford's labs at 
NCAR (Boulder, Colorado) and the results of the impactor data 
transmitted to us. Only results from the Yuma I and Yuma II flights 
will be presented here. 

Figs. ( 3 - 1 3 » a-h) are computer-produced size distribution 
curves for all measurements. Errors are indicated as dots on 
either side of the curve. The total number of particles measured 
in the . range of 0. 31 to 6.3 u.m was: 

Total 

Fl i ght A 1 1 i tude ( f t . ) Particles (Cm ) % Error 


Yuma 1 

1,000 

0.359 

6.8 


5,000 

0.125 

9.2 


10,000 

0.520 

6.7 

✓ 

40,000 

0.250 

6.9 

Yuma 1 1 

• 1,000 

1.03 

6.3 


5.000 

1 • 37 

5.4 


10,000 

0.252 

1.3 


40,000 

0.319 

6. 1 

Even 

though the particle count 

for Yuma 1 

1 is s i gn i f i cant 1 y 


higher than for Yuma I, both counts represent one effectively 
clear atmosphere. 

3.2.2 True Heading Problems: Relative Solar Bearings. - 

Since the polarization curves are defined with respect to the sun, 
it is necessary to determine the actual relative solar bearing, 
i.e., the difference between the true heading of the aircraft and 
solar azimuth. The solar azimuth at the relevant time, latitude and 
longitude, as well as solar elevation may be calculated using 
standard formulae and data from the A i r Almanac for the dates of 
interest. 

The video-tapes contain an accurate record of the aircraft 
heading relative to the sun as the instrument and TV camera 
scanned through the sun on the sky 1 i ght part of Leg 1 scans (1.25" 
on the screen - 1.5°) and oh the skylight scans during the low al- 
titude direct sampling runs. (The roll of the aircraft and solar 
elevation must be taken into account in translating the angular 
position of the sun measured on the TV screen to the relative 
solar bearing of the aircraft.) 
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Fig. (3-13) -- Direct Sampling Results. The following figures are 
computer produced plots of the number-size distributions of the 
col lected part iculates . (a) through (d) are from Y1 (12/29/71) at 
1,000, 5,000, 10,000, and 40,000 ft., respectively. (e) through 
(h) are from Y2 (12/30/71) at the same corresponding altitudes. 
The vertical coordinate is dN/d (log R) . See next page for (e) 
through (h) . 
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In correlating this information it was noted that there were 
discrepancies of as much as 6° in the true headings recorded by 
the. ADDAS. Fortunately* the navigator (R. Morrison) had made sex- 
tant sightings at the beginning and end of each leg. Comparing 
his records with the ADDAS printouts revealed discrepancies of an 
equivalent degree. 1 1' was decided to use the average of the nav- 
igator's sightings to designate the true heading for a leg. (there 
appears to be no consistency to the discrepancies) since the print- 
out values were clearly unreliable. 

In spite of the skepticism generated by this* we assumed that 
all other aircraft parameters on the ADDAS printout such as latitude, 
longitude, roll, and, pitch are reliable. The only other parameter 
which we could check was the pitch angle, which was averaged over 
the leg and used to calculate the scan angle. (Though the average 
pitch is recorded on the title blocks of the data curves presented 
later, it has already been incorporated in the scan angle, the ab- 
scissa of the graphs.) The elevation of the sun measured in this 
manner (only on Leg l's, of course) is always within a few tenths 
of a degree of the calculated solar elevation. This verifies not 
only the accuracy of the pitch, but the accuracy of the scan angle 
calibrations as well. 

Analysis of the solar passes, particularly on the direct 
sampling runs where the aircraft flew alternately into and away 
from the sun, indicates that the rotation axis of the scan platform 
is not normal to the true. heading axis of the aircraft, but canted 
at least one degree to the left (-1°). This has been incorporated 
in the relative solar bearings recorded on the title blocks of the 
data graph. 

3.2.3 Determination of Calibration Factors - Pre-f 1 i ght 
Elliptical Distortion . - Determination of the calibration factors 
prior to flight were made on the roof of MS using the rotatable 
col lar described in Section 3. 3 bf' the Appendix. The principal 
factors determined in these experiments are given in Table 3-3(a). 

The method is essent ial 1 y' ( b) of 3*1.1 with the instrument rotated 
through discrete intervals of 15°. 

In the course of these experiments, it' was discovered that 
instead of a circle, the data curve in the T - pl an e after taking 
f] and f£ into account was distorted elliptically with the major 
axis in the "2? = 135° (or 2 y = -45°) direction. Also, the angular 
intervals calculated for 2y varied by as much as ±2.2° from the 30° 
they should have been. As has been shown in Case II of 3.1*1* this 
type of distortion might be due to prism misalignment. Analysis of 
the data shows that the distortion can be accounted for by taking 
Y= 1 0 1 1 1 . Indeed, using a corresponding cal ibration matrix res- 
tored the curve to a circle and the angular intervals to their 
proper amount, and such equations were incorporated into the pro- 
gram. (These equations were later removed for reasons to be given.) 
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Attempts to make a direct measurement of the m i sor i entat i on 
were undertaken but the error in determining the or i entat ions 
was +0.5° for the optical bench setup used. It was thought that 
the measurements, though marginal, ruled out m i sor i entat ion of 
the prisms and the next source suspected was the electronic 
processor. Later tests described in Section 3*0 of the Appendix 
ruled out this possibility. (In preparation of this report it 
was discovered that the original calculation, indicating a value 
of 1°30' was necessary to explain the effect, was in error. The 
value for y given above restores mi sor ientat ion as a possibility.) 

Rinq-a-round at Moffett Field . - Originally, it was planned 
to use the true heading data from the ADDAS in assessing the ring- 
a-round data, but all other discrepancies aside, the updating 
proved inadequate for this purpose. Instead, the data was plotted 
in the Tl-plane and the center of gravity determined graphically, 
though again some arbitrariness was involved. Examples of the 
plots for C2 done with the computer are shown in Fig. (3-1*0. The 
values deduced from these graphs are given in Table 3, "3 (b), (c) 
and (d). The reliability of the values determined in this manner 
may be seen by comparing Table 3-3 (d) with (e) and (f) which 
were determined after returning to UCLA as described below. Note 
that there is at most 0.3% difference in the same f's and typic- 
ally 0.2% or less difference. 

Unfortunately, as will be seen we cannot believe that these 
values represent that accurately the value corresponding to the 
inflight data. This is indicated by comparing the RAR's at Moffett 
Field. Table 3-3 (b) values are fairly close to (d) values, but 
the values in (c) contain major differences, particularly in the . 
f] row and in f £ of C4. 

We might explain the differences in fj as being due to a 
large piece of dirt on the external window of the B barrel or on 
the image aperture lens or fibre optic bundle aperture, for B in 
the chopper compartment, while the singular departure of the f£ 
of Ck might be due to a piece of dirt (flake of paint on the 
appropriate sector of the C4 filter). 

In the course of these ,f-value determinations it was realized 
that the presence of the equations correcting the elliptical dis- 
tortion was fouling things up (as may be seen from Section 3- 1.1 
translation must be done first, then rotation or deformation) and 
the equations were removed. 

They were never replaced for two reasons: 

1) Due to the uncertainty as to the origin of the distortion 
it was felt we might be forcing a constant error to replace a 
random one. 
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2) If the distortion wasi caused by m i sor ientat ion of the 
prisms, it would have a minimal effect when the polarization was 
along one of the axes ( T)i or T)2 ) as is the case with the major 

plane of interest, the sun plane. 

In-flight Data Determination. - An attempt was made to obtain 
information on the calibration factors prevailing during the 
actual data taking by using method (c) of 3« 1 • 1 » the special 
angle method (SAM;. This method relies on the knowledge that the 
plane of polarization when viewing the zenith is normal to the 
sun plane. Analysis of the data on this basis yielded inconsist- 
ent results. 

Although the true heading discrepancies could account for' 
some of the inconsistency (a 1° error in heading can typically 
result in 1% error in the f-value), it is felt that the presence 
of dirt in the chopper compartment was the major source of the 
differences. 

Post-flight Studies . - On return to UCLA, the scan platform 
was mounted on a frame with casters. A pivot point, which could 
be lowered manually to contact the ground was located near the 
doorknob (instrument compartment) end. On the few occasions when 
sky conditions were suitable, the scan platform was taken down 
to the loading dock behind Si ichter Hall for further calibration 
factor studies. The pivot point was lowered and by slowly 
pushing the far (electronics box) end of the scan platform, a 
ring-a-round was performed. Table 3~3(e) shows that even though 
poor sky conditions prevailed for the Yuma II RAR, the factors 
determined graphically are consistent with those determined under 
far better conditions. (The area of the loading dock used was 
exceptionally flat, deviation from the zenith axis during rotation, 
excluding pebbles, being less than +0.2°.) 

The values in Table 3"3(f) were determined using SAM to 
study the validity of that method. The instrument was oriented 
at 0°, 45°, 90°, and 135° successively, using the edges of the 
scan platform as a sundial. This method yields knowledge of the 
plane of polarization to better than +0.5* Note that the values 
determined in this manner differ from those in Table 3 _ 3(e) by 
less than 0.2% on the average, indicating that this method would 
indeed be useful for determination from in-flight data if the 
true heading, or better, the relative solar bearing can be meas- 
ured to a sufficient degree of accuracy while the instrument views 
the zenith. This method should be kept in mind'for future flight 
experiments. 

On only three other occasions were sky conditions suitable 
for studies on the loading dock though on two days these con- 
ditions lasted only for an hour or so, forcing some hurried and 
curtailed measurements. No conclusive numbers were gained but 
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the following observations were made: 

1) What was initially thought to be a dramatic effect of 
the non-symmetr i c angular acceptance of the B and D barrels was 
traced to improper orientation of the instrument so that its axis 
was not aligned with the zenith. 

2) This misalignment came about due to the removal of the 
instrument package from the compartment and its reinstallation 
without checking the "level" of the instrument with reference 
levels on the electronics box. (There are no locating pins to 
insure accurate orientation on reinstalling.) The error in the 
direction of scan rotation could at least be the one-half degree 
observed. It was also noted that the axis of the instrument was 
off normal to ; the rotation axis by 0.2°. 

3) This uncertainty in orientation mitigated against any 
study of the maxing of Stokes "components due to multiple reflec- 
tions between the lens of the photopol ar imeter and the external 
quartz windows. Such an effect would be critically angle depen- 
dent. 

This discussion may seem somewhat lenqthy but has been 
deemed neces story to illustrate the difficulties in determining 
the calibration factors. 

For want of justification in using other values, those in 
Table 3-3(b) were used for reducing the data of Yuma I and those 
in (c) were used for reducing the data of Yuma II. It makes no 
sense to interpolate the f's since they represent ratios of 
transmissions. Once data plots were obtained, it was easy to 
see why inconsistencies had occurred in the attempts to obtain 
good values for the calibration factors. 

3.2.4 Experimental Results - Preliminary Discussion . - The 
reduced data for- - the rosettes of the Yuma I and Yuma II flights 
are presented in the following graphs (Fig. (3"15))« These graphs 
are' not assigned individual, numbers since they will be referred 
to by the code (see Table. 3-2) in the upper left hand corner, also 
given in the title block. (A complete set of graphs is presented 
only for C2 with selected examples for the other color channels. 
The reasons for this are given in the following remarks.) 

Several main features should be pointed out: 

1) The scan is from 30° below the nose of the plane through 
the nadir at 90° to 150° while tracking the target area. The 
scan is then moved at an arbitrary rate to the rear horizon (180°) 
and set at a uniform rate (about 2°/sec) to scan the skylight 
passing through the zenith at 270° and ending at the forward hor- 
izon at 360° (Fig. (3-8)). 
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Table 3.3 -- Calibration Factors 





Cl 

C2 

C3 

C4 

(a) 

Instrument alone 


0.9123 

0.9804 

0.9653 

0. 9701 


prior to f 1 ight. 

f 2 

0.9494 

0.9868 

0.9724 

0.9976 

* 

(b) 

Salt Lake RAR 

f l 

0.9200 

0.9825 

0.9675 

0.9750 



f 2 

0.9425 

0.9825 

0.9650 

0.9925 

(c) 

Yuma I RAR 

f l 

0.8970 

0.9735 

0.9565 

0 . 965 0 


1 . 

f 2 

0.9420 

0.9825 

0.9650 

1.0020 



r : 

' l 

n no nr 
0 . 

n nO C 0 
\j • yJ j 

0.9715 

■ 0.9775 

(d) 

Yuma II RAR 




f 2 

0.9488. 

0.9863 

0.9713 

0.9958 

(e) 

Loading Dock RAR 

f l 

0.9238 

0.9888 

0.9738 

0.9813 


2/15/72 

f 2 

0.9475 

0.9869 

0.9706 • 

0.9975 

(f) 

Load i ng Dock 

f l 

0.9282 

0.9918 

0.9740 

0.9839 


Special Angle, 
3/1/72 

f 2 

0.9468 

0 . 986 1 

0.9698 

0.9964 


(b) to (e) were with the instrument in the scan platform with 
external windows in place. 
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.2). For convenience in presentation as well as the conven- 
ience in data processing, the ground and skylight scans are not 
separated but presented in one continuous curve. 

3) the region of 160° to 1 80° represents the scan through 
the wing and tail section of the CV-990 and no significance can 
be attached to the peaks in that region. 

4) The region from 320° to 3^0° in LI curves is wherethe 
instrument scans through the sun. Such steep intensity gradients 
occur in the region of the sun that in a given quarter (color) 
cycle the 1 1 s relating to the outputs b„, b,, and b^ are all dif- 
ferent, rendering the values of polarization meaningless. 

5) Saturation occurs when passing through the sun, but such 
singular points are not marked on the graphs. The black stripes 
that appear on curves for C3 and C4 indicate where saturation 
occurs for extended periods, rendering the polarization values 
given meaningless. This saturation is an amplifier saturation 
(or equivalently exceeds the capacity of the A/D converter) not 

a 'detector saturation. The optional preamp gain was not utilized 
because, first, it was virtually impossible for the operator to 
know when saturation was occurring (only someone monitoring the 
visicoraer could tell when the system was saturating, and although 
spot-checks were made on the preliminary flights, there was in- 
sufficient personnel on the Y1 and Y2 flights to man the visicorder 
continuously) and second, the gain change would have sacrificed 
Cl and C2. 

6) Data taking on L4 of R2 1 s was invariably curtailed due 

to flight zone restrictions. . 

7) Dips and irregularities in the skylight portion of the 
Y2 curves indicate the presence of cirrus clouds on that day. 

; The major problem in the data indicated on" the curve by "A" 
will be referred to as Type "A" glych. Particularly strong examples 
are ‘found in Y2R2L1C2, Y1R2L3C2, and Y1R2L1C2. Not all such 
glyches are indicated and the reader will be able to identify many 
for himself. Type "A" glyches are characterized by a sudden 
break or transition from one continuous curve to another, indi- 
cating a sudden change in a calibration factor occurring on one 
color only, . In tracing these glyches in polarization back to 
the full data printout, it was found that they always correspond 
to a sudden.' change in one parameter (of bQ, bj, b^ ) only.. 

.• jf this glych were the result of a sudden change in the 
electronics, it would occur in all colors simultaneously, not in 
one co.lor alone. For the same reason, the sudden appearance of 
a hunk of dirt on the external windows is ruled out. The only 
plausible explanation for the behavior of this glych is the 

72 



movement of a relatively large piece of dirt, like a flake of paint, 
onto (or off of) the filter wheel. To verify this, the detector 
was removed and such a piece of dirt was discovered on one of the 
filters. The effect of such a hunk would depend on its size and 
position on the filter, but it is easily seen how the effect 
occurs in only one parameter if it is remembered that each barrel 
corresponds to one 120° sector of the filter. It was also verified 
that the dirt remains in place during many cycles of the instrument. 
While in flight, the vibration and changing orientation of the 
instrument during scanning could combine to shake it loose. 

At any hate, the behavior of the Type "A" glych indicates 
clearly the origin of difficulties in determining the proper 
calibration factors at any given moment and severely restricts 
the reliability of the data. We can, however, make a qualitative 
assessment of the reliability of a set of curves for a given color 
rosette, and altitude by observing how close the neutral points 
in the skylight scan come to zero (or less rigorously the corres- 
ponding point in the ground scan) on LI curves and how well the 
full set, i.e., LI, L2, L3> and L4, cross over at the nadir and 
zenith points. Since C2 fares better than the other colors in 
this respect and shows no saturation, it has been given the most 
attention. 

Before proceeding to the next section, it should be mentioned 
that the waves in the ground scans of the Y2 curves, particularly 
noticeable in Y2R2L2, are real ground effects. In reviewing the 
video tapes for this leg, dark, stripe-like diagonal patches 
drifted in and out of the photopol ar imeter 1 s field of view as the 
plane drifted right and left. Due to poor contrast and the 400 
Hz interference, it is impossible to discern the nature of the 
patches, i.e., whether they represented shadows due to slopes, 
vegetation, or some other coloration. 

Analysis of the Data . - The dashed curve in Y1R1L1C2 and 
Y2R1L1C2 is the theoretical internal field calculation for a 
homogeneous Rayleigh atmosphere. This calculation discussed in 
the theoretical program report section is for \ = 0.44mu , A = 

0.09 at 10,000 ft. altitude and the solar elevation (34°) or 
solar zenith corresponding to that of the experiment. The upward 
(ground) radiation and downward (skylight) radiation are presented 
in one continuous curve to conform with the data presentation, 
even though tne calculations were formally separated. 

The wavelengths for the computer calculations were chosen 
in order to ma!<;e use of existing tabulated information on re- 
flection and transmission matrices. The curve presented is 
interpolated from those calculations to the 0.44 mu wavelength. 

The albedo of 0.09 was chosen on the basis of the peaks of the 
data curves. For the given wavelength and type of ter ra i n*. ( dessert) 
the value is quite reasonable (ref. 30 ). 
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Fig. (3-15) -- Computer Plots of Polarization vs Scan Angle for 
Rosettes of Yuma I and Yuma II Flights. The following pages 
contain computer plots of selected data legs from the flight 
experiments. A complete set (16) of graphs for C2 is presented 
followed by one example for each of the remaining color channels. 
See text. (Section 2.4.1) for detailed explanation. 
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Note the excellent agreements particularly on the ground 
scan portions for both days showing that the aerosol has an in- 
significant effect on the polarization as would be expected for 
the particle counts' presented in Section 3-2,1. The departure of 
the theoretical curve on the left hand side of the peak in the 
skylight portion is not understood at present. It should be 
noted, however, that the skylight data from both Y1 and Y2 are 
in excellent agreement with each other except for the cirrus 
cloud dips. The agreement might be even better were it not for. 
the occurrence of a Type A glych at 200° in Y2R1L1C2. 

Though not shown here, curves calculated for the radiation 
emerging from the top of the atmosphere do not yield the same 
agreement, illustrating the importance of the internal field 
calculations. It was hoped that Dr. Takashima would be able tg 
calculate some curves for relative solar bearings other than 0°, 
but his hasty departure for Japan precluded this. Further calcu- 
lations will require both more time and money, both of which are 
exhausted at present. 

The entire set of curves from Yl and Y2 are presented here 
in the hope that workers either here or elsewhere will be able 
to make further analysis of them. 
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4.0 SUMMARY AND RECOMMENDATIONS FOR FURTHER WORK 


4.1 Instrumentation . --A detailed critique of the instrument 
is presented in the appendix and in Section 3*2.3 of this report. 

To summarize those comments* 'the instrument requires refurbishing 
before it can be used for further experiments* some of which is 
desirable and some of which is necessary. Once the sources of. 
systematic error are eliminated or reduced to an accountable state* 
the system will be capable of the high quality measurements 
desired. Evidence for this is seen in the low scatter of the 

data presented in the previous section, in which point-to-point 
connection results in a remarkably smooth curve. 

In many of the laboratory and field experiments the accuracy 
of the measurements was found to be limited by the three digit 
capacity of the A/D converter,. Since both the converter and 
tape recorder are presently outmoded and were the source of 
major problems, it is advisable to replace, these two components 
for future experimental work. 

The results of the flight experiment have demonstrated the 
validity of the instrumental system and techniques, and the 
feasibility of using a similar system for satellite measurements 
though some modifications may be necessary. A considerable 
practical advantage was gained in using the solid-state detector, 
thereby eliminating the need for a precision high-voltage photo- 
multiplier power supply. If, however, it is decided to is e a 
substantially decreased field of. view for the satellite instrument, 
it may be necessary to return to a photomultiplier as a detector. 

If the field of view is decreased by one order of magnitude or 
less, the solid state detector may be retained by operating in a 
biased mode. Though this would result in an increase in statis- 
tical errors, it would still be preferable to a photomultiplier. 

A common pitfall of experimental science is improper or in- 
adequate calibration of the instrumentation. Calibration prob- 
lems peculiar to the multi-barrel polarimeter as well as the 
general lack of polarimetry standards have been discussed. The 
need for exhaustive testing and calibration of the instrument 
prior to launch as well as methods for calibration checks in space 
(see, e.g. , S046FR (32)) cannot be overemphas i zed. 

4.2 Experimental Results . - For reasons already presented, 
the experiments were conducted under zero turbidity conditions. 

The comparison with the Rayleigh atmosphere calculation demon- 
strates the excellent agreement between theory and experiment 
that may be obtained with good instrumentation and good experi- 
mental procedures. As always, it is important to study first 
the simplest situations before haphazardly attacking the most 
complex situations. In this case, it means exhaustively studying 
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tracked homogeneous target areas such as the Gran Desierto 
primarily in the sun plane, under a variety of turbid conditions 
with simultaneous direct sampling. If this is not done, there is 
small hope of understanding results from a satellite where of 
necessity the "targets" may not be homogeneous (except of course 
over oceans, extensive cloud cover, or large desert areas) and 
the scanning plane may not be normal to the earth much less the 
sun-plane. 

As it was impossible to fly the instrument under turbid 
conditions, we have relied on the recent work of Rao and Takashima 
to demonstrate that a total atmospheric turbidity may be deter- 
mined from the polarization of emergent radiation. The tracking 
procedure and improved instrumentation developed under this 
contract should lead to greater accuracy and confidence in the 
turbidity determination. 

A. 3 Theoretical Advances . - The theoretical work conducted 
under the current program has produced two major advances in 
treating scattering by turbid atmospheres. 

The first and most general result is the equivalent charac- 
terization of the number-size distribution of polydispersions. 

This eliminates the multiplicity of distributions formerly used 
and reduces the number- of parameters needed to characterize hazes 
or aerosols. It is clear that this represents both a conceptual 
and’ computat i onal advantage in compiling a model library. 

The second significant result is the separation of the forward 
peak which will resul/t in a more tractable formulation of the 
computational problem leading to speed i errand therefore cheaper 
calculations. A corollary result is that the number-size distri- 
bution may be determined by study of the solar aureole from the 
ground. Though not directly applicable for the satellite experi- 
ment currently envisioned, this result may be used for ground 
based experiments in support of the satel 1 i te experiment. The 
result might also be used in a satellite experiment studying the 
aureole through the atmospheric limb to determine stratospheric 
particulates. Such an experiment does not require polarimetry 
but a polarimeter may be used if the field of view may be reduced 
to a size much smaller than the solar disc (~i° ). In this case, 
there would be no problem with signal strength since the intensi- 
ties involved would be many orders of magnitude greater than 
those involved in the earthward- 1 ook i rig polarimetry experiment. 

A. 4 Further Work . - The current state of the project suggests 
several avenues to be pursued in the near future. The outline 
below is somewhat arbitrary but is constructed to present those 
items which should receive the most immediate attention. Some of 
the items may be approached concurrently, others of necessity 
consecut i ve 1 y . 

So 



OUTLINE FOR FUTURE WORK 


Theoret i ca 1 

1) I ncorporat i on' of the theoretical results (separation 
of the forward peak) into model calculations. 

2) Formulation of plans for the model library including 
types of aerosol distribution and ground reflection 
characteristics as well as computational methods to 
be used. 

3) Compilation of tables and graphs to constitute a 
1 ibrary. 

Exper imental 

1) Further flight experiments on the Conva i r 990 with 
the photopol ar imeter over a particular area under 
various turbidity conditions with' concurrent direct 
sampl i ng. 

2) Ground-based aureole studies to determine number- 
size distribution. Since direct sampling would also 
be useful here, it would be advantageous to conduct 
such an experiment concurrent with the Convair 990 
experiment or with the other experiments in which 
direct sampling is being conducted. 

General 

1) Formulation of requirements for the satellite 
experiment including orbit, f ield-of-view, tracking 
accuracy, preferred target areas, etc. 

2) Development of the satellite photopol ar imeter and 
calibration methods. 

3) Formulation of plans for rapid reduction and assess- 
ment of satellite data (data management). 



4.5 CONCLUSION 


One of the original motivations for this study was the 
problem of whether or not increased industrial activity in recent 
years is causing an increase in particulate concentrations in 
the atmosphere on a global scale over and above that produced by 
naturally occurring processes. It is clear from the material 
presented here that it is entirely feasible to monitor global 
turbidity from a satellite photopol ar imeter on a long term basis 
to provide an answer to this question. 

: Application of the satellite data to yield information on 

the surface state of oceans, ground moisture, and crops, the . 
nature of cloud cover (ice or water clouds, etc.) as well as 
more detailed information on aerosols and ground reflection 
characteristics will depend upon further theoretical work and 
the availability of satellite data for analysis. 
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Test Reports on the UCLA/GE Photopol ar imeter 

/ 1.0 I nt roduct i on. - The important experimental objective of 

this project has been the development of instrumentation suitable 
for satellite installation (at least in principle) and capable of 
obtaining reliable data of sufficient accuracy. For this reason as 
well as for contractual ones, details of tests and analysis of the 
photopol ar imeter 1 s behavior are separated from the body of the re- 
port into this appendix. Some redundancy occurs and is desirable; 
it is hoped that this procedure will underscore rather than mini- 
mize the role of the PPM in the flight experiments and data inter- 
pretation. 

The instrument used in the current phase (i.e., contract) is 
based on ideas and prototypes originated at UCLA by T. A. Hariharan 
and W. L. Stevens under NASA grant No. NGR 05-007-041. ( ref . 28 ) 

The current design was developed in cooperation with the General 
Electric Space Science Laboratories at Valley Forge, Pennsylvania 
in breadboard form under an earlier subcontract (S046) and the model 
used in the flight experiments is a refinement of that model con- 
structed by GE under subcontract in the current contract. 

The final report on the S046 subcontract (S046FR) ref. 31) 
contains elaborate descriptions of the photopol ar imeter design 
which will not be repeated here. (The phenomenological design 
principles have been discussed in the body of this report.) A 
number of design changes were incorporated in the engineering 
model, however, the most significant being the i ncorporat i on of a 
solid state detector (United Detector Technology, Inc., PIN-10) 
instead of a photomultiplier. This resulted in an improvement in 
ruggedness, eliminated the need for a highly stable high voltage 
supply, and reduced space and weight requirements for both detector 
and power supply. 

Other significant differences are: 

1) A calcite compensating block in the I barrel with optic 
axis oriented similar to that in the G 1 an-Thompson prisms in the B 
and D barrels, i.e., normal to the optical axis of the barrel. 

This is either parallel or normal to the polarization transmission 
plane of the B prism. 

2) The fibre optics are not mixed at the detector or exit 
end of the bundle as in the S046FR studies. Each input occupies a 
120° segment of circular exit aperture. 

3) The gains of the AGC (Automatic Gain Control) are 1, 10, 
100 in the nomenclature of the UCLA workers where the signals are 
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called I, 2B, and 2D. The GE workers refer to the signals as 
2 , B, D and the gains as 2, 20, 200. In the reports by the GE 
workers transmitted here, their own nomenclature is preserved. It 

is hoped there will be no confusion on the reader's part. 

4) The G 1 an-Thompson prisms contain an interface layer of a 
new epoxy glue rather than the traditional Canadian Balsam. This 
was done to enable the instrument to endure lower environmental 
temperatures. 

5) An optional (lower) gain in the preamplifier was included. 

Section 2 of this appendix presents tests and results perfor- 
med at GE both prior to and after assembly, with the names of the 
responsible GE workers appended. For the sake of coherence, miscel- 
laneous reports have been edited and compiled. Section 3 includes 
comments by UCLA workers on tests observed at GE in May 1971 > (These 
are referred to as "acceptance tests" though in reality they were 
only demonstrations) as well as reports on pre-flight, post-fl ight, 
and in-flight behavior of the instrument. 

Before presenting these reports it is essential to point 
out some of the difficulties in evaluating photopol ar imeters in 
general, and this instrument in particular. First, there are no 
standard polarization sources. Textbook definitions of polarization 
parameters involve procedural prescriptions in terms of the meas- 
urement apparatus. Instruments such as this one may become the 
polarization standards of the future. Even sources for zero 
polarization do not exist; for some instruments the suri might be 
adequate, but only for those with fields of view less than half a 
degree. Devices’ for magnifying the angular field may well introduce 
residual polarization and could only be evaluated in terms of a 
reference photopol ar imeter or must satisfy some rigid theoretical 
conditions. Integrating spheres such as those used by both GE and 
UCLA described in detail in S046FR in Section 2.6 may be suitable if 
the observations are strictly on axis, but observations off center 
may detect a residual polarization. (A figure of merit for these 
spheres might be defined as the ratio of the area of the aperture 
to the hemispherical area, and degrees of polarization of this 
order can be observed off center. ) 

A unit polarization source can be defined by passing an arbi- 
trary source through a G 1 an-Thompson prims. (Note: Older definit- 

ions of polarization included use of Nicol prisms, but it is now 


a For convenience, the original Q, U, I, etc., notation used in 
S046FR is retained in this appendix instead of the more complex 
generalized notation intruduced earlier. 
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known that G 1 an-Thompsons are superior, the non-normal surfaces of 
Nicol permitting residual elliptical polarization.) For some pur- 
poses high quality polaroid sheets may be used, but the transmission 
of polaroid can vary significantly over the areas involved. 

This particular instrument has three independent and physically 
separated optical barrels. Even assuming perfect borcsighting and 
identical fields of view (3°) the fields of view do not overlap at 
laboratory distances. Thus, in tests involving integrating spheres, 
each 'barrel is viewing wel 1 separated portions of the target surface, 
and when a polaroid sheet is used, each barrel looks through a 
different portion of the sheet. Similarly, all three barrels 
cannot look through the same (external) polarizing prism simultan- 
eously. Rotating or translating the instrument to permit each 
barrel to look through the same prism or, portion of polaroid or 
portion of the integrating sphere target- is limited in accuracy by 
the accuracy with which the movements can be accomplished and by 
the stability of the light source. Also, laboratory sources have 
not provided sufficient intensity in the short wavelength channels 
and skylight observations have been necessary for studies of those 
channels. Certain evaluations have been made only for the 5000 A 
and/or 5800 a barrels and the results must be extrapolated to the 
others . 

Another significant problem in working with the instrument is 
that the system is not designed to give real -t ime readouts of data. 
Methods used for real-time readouts (quantitative) were limited by 
the equipment available as well as by peculiarities of the system 
itself and generally required measurements of one parameter at a time 
with significant time required to move the sampling pulse to another 
parameter (i.e., one Stokes parameter of one color, etc.). It 
was found that measurements recorded on tape and then reduced via 
the computer were more reliable and accurate, but this required 
excessive time delays , usually a minimum of three days. This 
presented severe limitations in working with the instrument. For 
example, if the skylight measurements were being made and the re- 
sults indicated either correction of some malfunction or alteration 
of experimental procedures, it was necessary to wait until the 
weather was again amenable for measurements. This was an extreme 
limitation in post-flight experiments. In the period from the 
beginning of January 1972 to the beginning of April it was possible 
to make skylight measurements on only two days. In order to begin 
preparation of this report, further experiments were abandoned. 

Not all of the tests performed either by GE or by UCLA are 
reported. In most cases, tests deleted were superceded by later 
tests or performance of the instrument in the flight experiments. 


*See Section 3 - 1 for a detailed description of such tests. 
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Even though the flight worthiness of the instrument has been dem- 
onstrated, the environmental tests performed at GE were an important 
part of the subcontract and essential for pre-flight tests. The 
letters and graphs verifying these tests are included as they were 
transmitted to us. 

The laboratory of the Atmospheric Radiation Group at UCLA 
does not have adequate facilities for performing of verifying 
many of the tests involved, : not only the environmental tests but 
those such as. cal ibrated sensitivity, angular response of the 
barrel s, , etc. Thus, the inclusion of the pertinent GE reports is 
essential . 



Table 1 -- Optical Component Tests 


Component F . 1 . 

Objective "B" 98 

Objective ”D” 100 

Objective ”1/2” 100 

Ap. Im. "B" 13 

Ap. Im. :: D !: 13 

Ap. Im. ”1/2” 13 


Prism ”B” 


Prism "D" 


Color-%T 


B i ref . ' 

1- 93 

2- 93 - 

3- 93 
A -96 

(2 . 9mm) 

6 < 10" 2 

1- 91 

2- 94 

3- 94 

4- 96 

( 2 . 9mm ) 

5 < 10" 2 

1- 91 

2- 94 

3- 94 

4- 96 

( 2 . 9mm) 

6 < !0' 2 

1 -96 

2- 97 

3- 95 

4- 95 



i -96 

2- 97 

3- 95 

4- 95 



1-96 

2-96 

3- 95 

4- 95 




FOV ( 63 28^) Re j ect . 


1-35 

16.3° 

< 10 5 

2- 42 

3- 42 

4- 46 

(16.7°) 


1-36 

15.3° 

< 10 5 

2-43 

(15.2°) 


3- 43 

4- 46 

' 
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Table 1 -- (cont.) 


Component 

F. 1 . 

Color-%T F0V(6328X) -Reject. 

Compensator “1/2” 


1- 73 

2- 81 
3-90 



4-94 

Fiber Bundles 


1-66 

2-66 

3-66 

4-63 



i 




i 
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2.0 Tests at General Electric : 2.1 Preassembly Optical 

Component Tests (T. Wise) . - Table 1 is a summary of the various 
parameters of the components. The leftmost column groups the 
components: Objective lenses, aperture imaging lenses, calcite 

prisms and. compensator block, and the fiber bundle. Following this 
there are tabulated the transmissions of the components at the 
four different wavelengths. The "Dispers" column lists the chrom- 
atic and polarized dispersion, in terms of change of focal length, 
for the lenses and compensator block, respectively. Under "Biref" 
there is listed a maximum value of birefringence for the objective 
lenses. This is given in terms of the induced optical phase shift 
b; the following formulas relate 6 to the residual birefringence 
(An) of a component and to the light intensity transmitted by the 
second of two crossed analyzers. 

I = 4A 2 sin 2 0 cos 2 0 sin 2 6/2 


Where 0 is the angle between optic axis of material and the axis 
of the first polarizer, A is the amplitude of the incident radi- 
ation, and 6 = 2jX. y dnn is the phase difference introduced by a 
birefringent material of thickness d^ This intensity will be a 
maximum for 0 = 45° where I = n sin z a/2. 6 less than 10“^ repre- 
sents induced spurious polarization of less than one part in 1 Cr 
at P = 5%- The "Fov" column in the data sheet gives the measured 
polarized fields of view of the two calcite prisms. The values 
given, however, relate only to the normal input face of the prism. 

As we learned later, the angles should have been referenced to the 
optical axis rather than to the mechanical axis. The prisms have 
a half degree deviation which makes the actual polarized field of 
view larger than this data would indicate. • The column labeled 
"Reject" gives a maximum value for a passage of the ordinary ray 
in the calcite prisms. My tests indicated that the prisms were 
at least as good as specified in this respect, possibly better. 

The basic set up of the test bench is as follows: An 11-inch 
refracting collimator was used to check the collimation of the light 
coming from a 3 inch f/1 mirror. All of the component testing was 
done with light from a gra i n-of -wheat lamp at the focal point of 
this mirror. Ahead of the lamp was an aperture which restricted 
the light bundle to the diameter of the component under test and 
ahead of this was a filter holder. I was able to get sufficient 
illumination at 3800 A by consuming several lamps at 20% over volt- 
age. To determine the dispersion of the lenses and of the compen- 
sator block, the same basic set up was used with a needle point 
and a microscope. A UV image converter was used for the 3800 A 
test. Placing the needle at the focal point of the lamp filament 
image removed chromatic effects in the microscope. ' 
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■ To check the component transmissions in addition to the mirror 
and the light source, a remotely operated aperture and shutter were 
used which permitted easy subtraction of noise present at low signal 
levels and hence provided accurate (+1%) values of light trans- 
mission. Readout was via a Schottky detector, opamp, digital 
voltmeter combination. A HeNe laser whose output was fed through 
the two G 1 an-Thompson prisms and onto a detector was used to de- 
termine the prism rejection and polarized field of view. The 
dynamic range of the signal levels went from 0.1 millivolt to 10 
volts, thus permitting resolution of 1 part in 1 0 5 when the uncrossed 
polarizer signal level was set at 10 volts. 

Figs. 1 - 4 are the filter transmission curves provided by 
Baird-Atomic. In the instrument, the 5800 A filter has been modified 
by the insertion of a No.. 8 Wratten filter to eliminate the UV 
wing in the Baird-Atomic filter. 
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Transmittance of Interference Filter for- Color Channel One (Cl) Nominally 




















Transmittance of Interference F i 1 ter for Color Channel Four (C4) Nominally 










2 . 2 Field of View: Bores i qht i nq and Angular Response of the 

B> D> and ? Barrels . - On June 7 , (1 97 1 ) the instrument was set 
up in front of the 100" collimator with a pinpoint of light as a 
probe to determine the boresight, cut-off. and uniformity charac- 
teristics of the instrument's field of view. The results of the 
first test, when graphed, indicated that the individual fields of 
view were within a tenth of a degree of the specified three degrees 
and were bores ighted with each other to within a tenth of a degree. 

Also, cut-off was as expected at the 5000 a wavelength (Fig. 5)* 

However, two oddities appeared in the B channel in particular. 

First of all, its overall signal level was 20 to 30% higher than 
either of the other two barrels. Secondly, it had a notable tail- 
off of the signal on one side. 

In an effort to understand what these effects were, the test 
was repeated with the i barrel at the 135° position, which has the 
effect of orienting the prism in the D barrel precisely as the 
prism in the B barrel had been oriented in the' previous test. 

Although this data appears in the notes, it was not graphed since 
we found that one of the difficulties in the measurement of the B 
barrel was a 10% residual polarization effect in the collimator 
favoring a greater signal for that channel. We concluded that the 
data as originally taken was stili valid for the determination of 
boresight, cut-off, and size of the fields of view. However, un- 
iformity was in doubt. 

With the purpose of checking the uniformity o-f field alone, 
the integrating sphere was set up 30 feet from the instrument and 
closed down to a 3 inch opening. The instrument was then scanned 
across this target with the I barrel at the 180°, 135°, 90°, and 
45° positions. This data appears in graphical form (Fig. 6) and 
indicates that a small light source does indeed produce a larger 
signal in the B channel; however, the difference now is only 10%. 

This 10% difference may be accounted for by the slightly restric- 
ted field of view of the B barrel, causing the downward trimming 
of the other two barrels in order to match integrated responsivity 
of the three channels. In addition, the tail-off of the signal 
in the B channel appears as much as it did in the previous test. 

In fact, there is a similar, but lesser effect in the D barrel 
when the D prism is oriented as the B prism was oriented in the 
previous test. 

The physical significance of this effect may be stated as a 
character i st i c of the G 1 an-Thompson prism and its small changes 
in refractive index across the interface. Internal reflection of 
the extraordinary ray increases as its incidence angle on the inter- 
face increases. This explanation has been confirmed by the actual . 
orientation of the prisms in each of the barrels and by the direction of 



scan of the point light sources in these tests. The operational 
significance of this is the greater contribution to the overall 
signal that a target on one side of . the field of view makes rela- 
tive to a target on the other side of the field of view. Since 
proper operation of the instrument demands uniform targets, this 
should be no problem at all. 
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Fig. 6 -- Results of Angular Response Studies Using 
rating Sphere with 3" Opening 30 feet away with the 
at (a) 180°; (b) 135°; (c) 90°; and (d) 45°. 


the Integ- 
I Barrel Set 
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2 . 3 Sensitivity: Linearity and Intensity Calibration of the 

I /2 Barrel ♦ - The sensitivity tests which were carried out early in 
June determine the system responsivity at each wavelength, the 
linearity over one decade in any one color, and linearity over 
several decades among all colors. A pleasing aspect of the respon- 
sibility-information is that the system at 3800 A performs three 
times as well as originally predicted by the minimum system design. 
This has the effect of contracting the overal 1 dynamic range of the 
output signals, since the signal at 5800 fi is no more than predicted. 

Fig. 7 shows the results of these tests with the output of 
the photopol ar imeter in volts plotted vs. the reading in amperes 
of the calibrated E, G, & G radiometer ( E, G, & G High Sensitivity 
Detector Head, Model #585-66 S/N 78). Table 2 gives the resultant 
calibration of the 1/2 readings taking into count the calibration 
of the radiometer. 
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2 . 4 Environmental Tests (L. Napaluch and T. Wise) . - On May 
19 and 21, 1971 the i ns t rument and electronics box were subjected 

to temperature and altitude, and vibration tests certified by the 
letters and graphs immediately following. During both tests, a 
light source was arranged and output from the processor was moni- 
tored (see Fig. 8). There was. considerable ripple in the lamp and 
its intensity varied as various units in the vicinity cycled on 
and off. In addition, the data was very noisy due to disturbances 
from the refrigeration unit heaters and the low altitude and high 
altitude pumps. Because the instrumentation must be grounded, 
neither el ect ron i cs . box nor optical assembly were insulated from 
the chamber, and a ground loop was established through the instru- 
mentation ground. For these reasons, noth i ng concl us i ve can be 
drawn from the data observed during the environmental test except 
that no inexplicable or grossly aberrent behavior was observed, i.e., 
the instrument and electronics behaved well -during the test. No 
difficulties were encountered in stopping and restarting the motor 
at 0°C. . 

Measurements taken with the integrating sphere before and 
after the environmental test were compared, the E, G, and G radio- 
meter being used to set up identical intensities in the sphere. 

The intensity signals changed very little; in fact, there is a 
uniform half percent increase in the post-env i ronmental test 
studies. This is within the prev i obs I y determ i ned . i ong term 
stability of the instrument. In many cases there were no discrep- 
ancies for the Q and U readings and where there were they represent 
maximum changes of 25% of the signal voltage. This is very good 
considering the fine balance that exists in the Q and U outputs. 
Moderate changes in the instrument and in the testing of the 
instrument will produce voltage changes on the order of decades. 

Subsequent to ( 1 ) specified vibration levels up to 5 G's and 
residual levels up to 9 G's (2) temperature soaks at 0°C and 50°C, 
and (3) vacuum levels down to j 00 torr, there was no detectable 
mechanical instability greater than 0.5%. 
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Fig. 8 -- Photograph Showing a Typical Set Up of the 
Photopolar imeter During the Vibration Tests at General 
Electric, Valley Forge Space Center. 
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2 . 5 Sources of Extraneous Noise (L. Napaluch) . - ( F rom May 
26 to May 28, 1 9 7 1 Mo Kopa aricFDr. Bradbury witnessed acceptance 
tests of the instrument at the General Electric Space Sciences 
Laboratory.) During the data gathering phase of the acceptance 
tests a considerable amount of noise or signal output variation 
was evidenced with the instrument looking into the integrating spher 
at essentially low polarization light. This variation was greatly 
in excess of that expected from the known pre-amp noise, so further 
exploration was indicated. 


The 

1 ) 

2 ) 

3) 

4) 

5) 

6 ) 


noise can be broken 
Ripple on the 1 i ght 
Pre-amp noise 
Pre-amp input cable 
Photodiode noise 
Processor drift 
False triqqerinq of 


down into the following list: 
source in the integrating sphere 

noise (microdot) 


the data system 


I tern 1 . - The light source ripple arises because the excitation 
is a full wave rectified output from an A.C. regulator. The light 
intensity ripple of 120 Hz is less than 1% and is somewhat atten- 
uated by the 3.6 millisecond time constant of the photopol ar imeter 
leaving a residual of about 0.2%. 

This source of disturbance is inherent in the instrumentation 
and makes it difficult to read the processor output accurately. 

To eliminate it, a new regulator, battery powered, has been built 
and tested. All further instrument tests with the integrating 
sphere will use this new regulator and hence will be free from 
ripple. 

Item 2 . - The pre-amp noise was the same as it had been when 
measured and calculated previously. This constitutes a limitation 
in photopol ar i meter performance at the highest gain, the gain with 
an open circuit link. This gain provides full output of 10 volts 
with 5*2 X 10 ' amps _4nput. When the pre-amplifier gain is set for 
10 volts at 5-2 X 10 'amps, the design point for the instrument, 
this noise will be only 10% of what it is at the highest gain. At 
the nominal design gain, the pre-amp noise will not be significant 
in determining instrument performance. 

I tern 3* - With the cable alone placed on the pre-amp input, 
the re “may or may not be an increase in noise, depending on the 
level of vibration the cable sees. When cable and pre-amp were on 
a table near the masonry wall, there was no significant increase 


in noise. When they were placed on the test table near the center 
of the floor, a place of known high vibrational level, the noise 
went up by a factor of two at the pre-amp output. As noted in 
Item 2, this would contribute to a limitation in photopol ar imeter 
performance at the highest gain. Likewise, noise contribution 
of the cable at the nominal design point would not be significant. 

I tern 4 . - The photocell contributed sufficient noise to 
triple the noise output of the pre-amp when connected through the 
microdot cable. Further, during the high temperature test a 
rather large increase of about 8 to 10 times was seen in the noise 
at +50°C. The photodiode was probably contributing most of this 
increase but this has not yet been established. Although at the 
nominal design gain this noise would be much reduced, it would 
be significant in limiting photopol ar imeter performance. 

I tern 5 . - From the specifications of the processor amplifiers, 
an expected drift is about +0.5 millivolts at the input. When 
the barrels were closed and this residual voltage (an intercept 
error and not a slope error) was observed, it was consistently 
below this level. However, under varying temperature conditions 
such as initial turn-on stabilization, it can exceed 0.5 milli-volts. 
The use of 40K amplifiers in place of 40J units (Analog Devices, 
inc.) helped keep this error to less than 0.5 millivolts as measured. 

This drift error can be a limitation in photopol ar imeter 
performance particularly at the nominal design gain. At the high 
gain pre-amp setting it contributes about equally with the noises. 

I tern 6 . - The open breadboard of the sample-hold circuit and 
its driver used to obtain data at a point in the processor output 
was subject to false triggering whenever certain equipment was 
operating in the vicinity. This caused a problem only when trying 
to average mentally many successive readings. 

3 . 0 Tests at UCLA - 3.1 Comments on the Acceptance Tests * 
at GE (R. Bradbury) . - The source for the tests observed was a 
quartz- i od i ne tungsten lamp mounted in a four-foot integrating 
sphere. The source as actually observed by the instrument was a 
5 or 6-inch diameter plane target painted with the same material 
as was the interior surface of the sphere. 

Many of the problems involved in evaluating the data are 
related to the source, e.g., for any one measurement, each barrel 
of the photopol ar imeter is looking at a different area on the 
target. Additionally, there are short term temporal variations 
of intensity of about 0.08% and long term variations within +0.5%. 

The short term fluctuations do not appear in the G.E. data 
records, since T. Wise recorded the average of a number of suc- 
cessive readings on the digital voltmeter. My own observations 
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of these fluctuations are summarized in the following typical 
data, taken with 0% polarization, 0° rotation of the photopolari- 


meter: 




Color 

Fl uctuat i on 

Gain 

Parameter 

5800 8 

+0.003 v 

2 

for l(4.29v) 


+0 • 6 

200 

for Q and U 

5000 % 

+0.002 v 

2 

1(1.360 v) 


+0.3 

200 

Q and U 

4400 A 

+0.020 

20 

1(2.725 v) 


+0. 20 

200 

Q and U 

3800 A 

+0.09 

200 

1(3-04 v) 


+ 0. 10 

200 

0 and U 

The 

to peak 
too 1 ow 

fluctuations in the 3800 $ signals represent 2 mv peak 
noise which is signal independent, since the intensity 
in this color channel to see the effect of the source 


f 1 uctuat ions. This noise, illustrated in Fig. 9, which shows 
oscilloscope traces of ten consecutive cycles of the Dhotopolari- 
meter, currently represents the limiting error source, far above 
the expected 0.5 mv error due to the processor electronics. The 
sources of the noise have been discussed further by Napaluch in 
Section 2.5* It is, however, small compared to the source ripple 
which in turn is small compared to discrepancies due to the spatial 
inhomogeneity of the source. 

The signals in both the 4400 $ and 3800 % channel are too 
weak to make any meaningful assessment of the operation of the 
instrument at those wavelengths, so my main concern in this report 
will be the 5000 A and 5800 A channels. Another source, possibly 
skylight, must be used to evaluate the short wavelength channels. 

Fig. 10 shows five consecutive traces of the output of the 
processor with the signals identified above. (The signal sequence 
is the same in all similar photos to follow.) The sampling pulse 
of the test circuit rather than the gain circuit is shown. The 
signal at this point (I in the 4400 A channel) is read out of the 
digital voltmeter. The situation is the same as for the table 
above. 

The effect due to the source ripple may be seen clearly in 
the Q and U signals in the 5800 A channel, where in both cases at 
lease four different values for both signals may be discerned. 


*This section is excerpted from a report on the acceptance tests 
prepared for Professor Z. Sekera. 


The effect is not observable in the corresponding I signal since 
the gain there is a factor of 100 lower. (Note: The gain factors 
in data taken by G.E. are given as 2, 20, and 200. In my own 
reduction of the data I read these as 1 , 10, and 100 to avoid 
the absurdity of dividing every piece of data by 2.) 

The effect due to the target inhomogeneity will be discussed 
in the following remarks. 

0% Polarization Test . - The laboratory set-up for the photo- 
polarimeter performance tests are shown in Fig. 11. The photo- 
polarimeter was mounted in a rotatable mount, not indicated in the 
figure. For lack of evidence to the contrary, we will assume 
that the axis of the mount is coincident with that of the photo- 
polarimeter so that, as the instrument is rotated, the area viewed 
by the i barrel, for example, generates a doughnut-shaped area on 
the target. 0° rotation of the photopol a r imeter indicates that 
the j barrel is at the topmost position in the rotatable mount. 

The B barrel is then viewing the 120° position on the target, i.e., 
the area viewed by the i barrel when the instrument is rotated 
through 120°, (this defings the direction of rotation) and the D 
barrel is viewing the 240° position. Thus the instrument 
notated by at 0°", "B at 120°", and "D at 240°" are equivalent 
statements. 

The trim factors f and g are defined on pages 124-125 in S046FR 
and a method for determining these factors given on pages 126-127* 
(See Section 3*1.1 of the report body.) Ideally, these factors are 
unity and should be the same for all color channels. Unfortunately, 
the method given in S046FR is unsuitable when using the integrating 
sphere as a source. The data from the no-pol ar izer^rotat ion 
tests are shown in Fig. 12 for the 5000 A and 5800 qA channels. 

The trim screw adjustments were made for the 5800 A channel with o 
the instrument at 0°. This corresponds to the 0° point for 5800 A 
being closest to the center of the graph. (Distance from the 
center is the degree of polarization.) The circles about each 
data point indicate the error limit due to the fluctuations 
discussed in the previous section. It is clear that the deviation 
of the points is much greater than could be accounted for by 
these fluctuations. Note that except for the 240° point in both 
cases the data points roughly describe an ellipse as the instrumet 
rotated. From page 127 of S046FR such an ellipse would be expected 
if the source were uniform, but with some small residual polari- 
zation. The center of the ellipse would be given by 1-f*, 1-g and 
the ellipticity would also be related to the trim factors. In that 


*At the time this report was prepared, the 90° mi sor i entat i on of 
the B and D prisms had not been discovered. 
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case* however, the point would travel around the ellipse twice as 
the instrument is rotated through 360°, so residual polarization 
could be only a small contribution to the deviations. This is 
confirmed also by the fact that the ellipticity does not corres- 
pond to the ratio of g/f, i.e., the figures would be more like 
circles to correspond to the location of the centers. 

As indicated in Fig. 12, it was found that it was impossible 
to trim the instrument well for all the color channels simul- 
taneously. This has been explained by the difference in thq 
transmission properties of the compensator calcite block in the 
i barrel compared to the Gian Thompson prisms in the B and D 
barrels (Table 1). Thus the 5000 A data is considerably displaced 
from the origin when the instrument is trimmed for 5800 A. 

Such an effect should be expected since the ratios of the 
trim factors for the different color channels are related to the 
barrel transmissions by equations such as: 

f ( X, ) T( X, » B ) . T(X 2 ,B) 

f(Xl) T( X, * I ) * T(X 2 ,I) 


where f(x ) is the trim factor f for color \ and T(x»X ) is the 
overall transmission of barrel X at color x> etc. 

With data like this, the best method of determining the trim 
factors is to rough in an ellipse to the data, and graphically 
determine the center. (This was done and the results used for 
reducing the data in the subsequent test.) 

Let us consider that the discrepancies are due solely to the 
snatial i nhomoqene i ty of the target. Assuming that the polarizat- 
ion is exactly zero and f and g are both unity, then it is 
read i 1 y seen that 


g(9) = 1(6) - 1(9 + 120°) 

i(6) 1(0) 

u( 6) 1(9 ) - 1( 6 - [20^1 

i(e) TT97~ 
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(Note: Following S046FR q(0) and u (0) designate the quantities 

as read out directly from the instrument when rotated by 0. 1(0) 

designates the intensity of the target at position 0 on the targeti 
i.e., as measured by the 5 barrel at 0.) 

Thu^ deviation from zero by .008 (.8%) for P would correspond 
to the same order of inhomogeneity. This could easily account for 
the distribution of points observed and conceivably the readings 
of the i barrel as the instrument is rotated might be taken into 
account in reducing the data. Unfortunately, the time between 
readings is too long and the long term source fluctuation leads 
to changes in the intensity of a comparable magnitude (+0.5%)* 

100% Polarization . - For the 100% polarization test, a sheet 
of HN-22 Polaroid film is taped to the front of the integrating 
sphere with its transmission axis horizontal, i.e., effectively 
parallel to the transmission axis of the B barrel. 

Fig. 13 shows a plot of the raw data from this test for 
5800 A and 5000 8, together with a circle for ideal 100% polar- 
ization. In this case, the gross deviations from the unit circle 
can be reduced by taking into account the source inhomogeneity. 

This works because the i nhornogene i ty is due not so much to the 
sphere properties as to the non-uniform transmission character- 
istics of the polaroid. On the basis of the data, this may vary 
as much as 6% around the doughnut. Fig. 14 with the data correc- 
ted according to the formulas 


1(e) = i(0). 


m 


. i(9) - q(0.) 
f i ( © +120) 



i( 0) - u( 6) 

g i ( 0 - 120) 
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shows considerable improvement. (The f and g values used are Q 
those determined from the zero polarization data.) The 5800 A 
data have an average . va 1 ue of P and 0.992 with RMS deviation 
0.0027 while the 5000 A data average is 0.985 with RMS deviation 
0.005. The relative deviation between the two color channels 
is compatible with the relative fluctuations (compare the circles 
about the data points in Fig. 12). The absolute error is in turn 
compatible with the ellipse of the no-pol ar izat ion data ( AP~+0. 007 ) ), 
indicating that this is the best that can be done in terms of the 
data presently available. 

One more comment must be added. It was impossible to obtain 
a ratio of less than/^10 ^ when measuring the extinction ratio of 
the B and D barrels using the polaroid sheet. Measuring the 
Polaroid's extinction ratio by crossing two pieces of film in 
front of the i barrel verified that this was due to the polaroid. 

The preceding discussion clearly indicates that the tests 
performed with the integrating sphere are inadequate to demonstrate 
the ultimate capabilities of the instrument. 
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Fig. 9 -- Ten Consecutive Traces of the Processor Out- 
put with All Barrels Blocked Off. The time from one end 
to the other is one second, the duration of the instrument 
cycle (see Fig. 10). This illustrates the limiting 
noise of the system. The scale is 0.5 v/cm and the 
instrument gain is 200. 
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Fig. 10 -- Five Consecutive Traces of the Processor Output 
with All Barrels Open and the Instrument at 0° Looking into 
the Integrating Sphere (no polaroid). The information 
sequence is shown above the photo. The lower trace shows 
the sampling pulse of the test circuit sampling here I of 
the 4400 A channel. The scope scale is 2 v/cm and the 
instrument gain factors for the signals are (reading from 
left to right) 200, 200, 200, 200, 200, 20, 200, 200, 2, 
200, 200, 2. 
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Block Diagram of Laboratory Set-Up for Performance Tests of the UCLA PPM. 









f ig. 12 Uncorrected Data From 0% Polarization Tests for 5800 
and 5000 A. 
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3.2 Power Supply Problems . - When the instrument was first 
delivered to UCLA by Dr. Frost a variable voltage power supply 
was used to test the instrument to insure that it had survived 
its transit. Even though great care was taken to prevent an 
overvoltage, the instrument failed to function the next time. It 
was determined that every transistor in the power converter 
section of the primary electronics box was gone, due perhaps to an 
overvoltage spike from the variable voltage supply when it was 
turned off. 

Construction of the proper external supply was rapidly 
completed and no further difficulties of this type were encoun- 
tered. A power failure in the course of the flights was corrected 
by giving positive retention to the one circuit board in the power 
converter section which did not have it. 

These comments are included not so much as a criticism of 
the instrument, but as a word of caution to its user. 

3.3 Integrating Sphere Tests : Linearity, Gain Accuracy , 

"Polarization 11 of I.- In order to have a source usable for rudimen- 
tary laboratory tests at UCLA an 18" integrating sphere with a 6" 
circular target suspended in the center by wires was constructed 
wilfi a quartz- i od i ne lamp as a source. A variable voltage regu- 
lated D.C. power supply was constructed for the lamp. No noise 
due to ripple wa s present. 

o 

Linearity of the I barrel for color channels 3 and 4 (5000 A 
and 5800 A, respectively) was checked using a set of calibrated 
neutral density filters in conjunction with the sphere. Some 
discrepancies were noted, but by recalibrating the filter with a 
Cary Spectrophotometer Model 14 it wa s verified that these were 
due to errors inthe calibration curves. These filters were also 
used to check the accuracy of the gain settings in the AGC. 

Within the accuracy of the filter calibration the gains are what 
they should be (10/1). 

The rotatable mount used by GE was borrowed temporarily until 
our own could be constructed. It was noted that due to a displace- 
ment of the optical axis (the "center of gravity" of the three 
barrels) from the mechanical axis (the center of the mounting bolt 
circle) the three barrels di^j not view identical areas of the 
sphere target under the +120° permutations discussed in Section 3. 

To remedy this, an eccentric collar was provided in UCLA's rotat- 
able mount. The mount was also provided with locating holes every 
15° and a spring loaded locating pin. The instrument could thereby 
be rotated quickly and accurately through integral multiples of 
15° as desired (see Fig. 15). This collar was used primarily in the 
trim calibration factor determinations by observing zenity skylight 
from the roof of the Mathematical Sciences Building (MS) at UCLA. 
This is discussed in detail later on. The integrating sphere is 
inadequate for determining the calibration factors, but once 
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these factors are known the instrument can be used to evaluate the 
sphere. In this manner, the polarization of the target region 
viewed by 1 at 0° was determined to be .05 + *01. (The primary 
source of the error is due to spatial variations in the target 
i ntens i ty . ) 

There had been some question as to whether the presence of 
the compensator block and the fibre optics might engender a 
polarization sensitivity in the I barrel. With the instrument 
viewing the integrating sphere at 0 , a polaroid was rotated in 
front of I, care being taken to insure that the center of rotation 
was on the axis of the I barrel. An intensity variation equivalent 
to .055 +002 polarization was observed for color 3. Thus, for the 
I barrel at color 3 the difference in transmission for the ortho- 
gonal polarization components could at most be 1%. It is felt 
that the accuracy of this test could be improved if time did not 
prohibit further testing. 

3.4 Zero Levels and Noise . - During the initial testing of 
the instrument on the roof of the MS building exceptionally 
large abnormalities were observed which were traced to pick-up of 
60 cps signals being radiated by some other roof-top experiment. 

It was discovered that the detector cable shield was ungrounded 
(according to design) and the situation was remedied by putting 
the microdot cable through a grounded copper braid shield. This 
solved most of the noise problems previously experienced. 

In subsequent tests of the instrument with all barrels 
blocked, it was found that while the RMS deviations were typic- 
ally 0.6 mv, the mean reading of a given parameter could be non- 
zero by as much as lOmv, the I's usually being highest. These 
zero levels (Z's) could change spontaneously. it was verified that 
the Z's could not be due to light leakage or scattered light but 
might possibly be due to static charges on the chopper wheel or 
filters. In most cases these Z's would lead to insignificant 
errors but they have been taken into account in data reduction as 
a matter of course. During data runs the barrels would be 
blocked before and after and the Z's averaged. This practice 
was followed both in ground and flight measurements, primarily 
as a check on the instrument behavior. It should be noted that 
the noise parameter given above is effectively the same as the 
0.5 mv error limit specified for the electronics. This excellent 
behavior was also observed in the flight experiment data. 

3.5 Comparison wi th Another Polarimeter . - In mid- October 
the photopolarinieter and a rotating prism instrument (Dr. N . C. R. 
Rao's) were run simultaneously for ten minutes while viewing the 
zenity from the roof of the MS building. The instruments were 
separated both vertically and horizontally by distances of about 
15 feet. At no time during this period did results from these 
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Fig. 15 -- Photo of Laboratory Set Up Showing the Photopolar imeter Mounted in 
Rotatable Ring Viewing the 18" Integrating Sphere. 




instruments differ by more than 3% of the polarization, the average 
difference being 1.8% for the 5000 A channel which was the only 
common filter region for the two instruments. 

3.6 Residual Charge Effects . - Some concern was expressed 

about the possibility of parameter mixing due to residual charge 
effects in the processing capacitor. Originally it was suggested 
that residual voltages from the usually high signal color 4 
channel might be carried over into the usually low signal color 
1 channel (The intensity difference is more exaggerated with 
laboratory light sources than with skylight radiations ) and 
a model processor was built to study such effects. It was real- 
ized that if such residual voltages might be carried over from 
one color channel to another, then the effect would even be larger 
within the parameters of one color channel where the time differ- 
ences were much shorter. 


For example, if after the Q cycle a fraction sQ remains on 
the processor capacitor, then 2 the U reading will be U + eQ and the 
I reading would be I + eU + e^Q^I + eU. Correspondingly, we would 
have 


71, = Q 11, - ell, 7], 

i + eu 


' _ U + cQ TL + eT| 

\ ~ I + eU l 


2 



I 


or, assuming we can ignore terms of "third order" (i). and TU 
are less than unity by definition) 1 1 

nj - H 

I 

ii2 = ^2 e i 
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This leads to an elliptical distortion indistinguishable from 
that due to prism m i sor i entat i on (See Section 1.1.3 in the body of 
the report). 

To search for such an effect in all color channels, a test 
was run viewing zenith skylight. When the B and D barrels are 
blocked the readings of the Q, U, and 1 outputs should be -I, 

-1, I, respectively. The results of this test showed clearly 
that within the usual fluctuations, the readings were identical 
(except for sign, of course) indicating no intrachannel parameter 
mixing. Actually, the elliptical distortion from such an effect 
would be in the opposite direction from that observed (See section 
2.3.1 in the body of the report). 

3.7 Black Paint: A Major Problem . - In late October the 

detector housing and detector were removed in order to investigate 
the sudden appearance of abnormally high Z's. It had also been 
noted that the detector was loose in the housing and could 
easily be turned in the process of connecting the detector cable 
(This can be avoided by disconnecting the Microdot coupling and 
not the BNC-to-Microdot adaptor) and it was absolutely necessary 
to retighten it. 

When removed it was discovered that the detector housing had 
been sprayed with black paint which was deteriorating and peeling. 
It was also noted that the edge of the housing had been grazed 
by the chopper wheel. The black paint was removed and a black 
oxide coating was plated on the brass housing. 

Before replacing the detector, it was necessary to separate 
the optical and mechanical modules to clean the remains of the 
peeled paint from the chopper compartment. Instructions for this 
disassembly were relayed by phone by Mr. T. Wise and the use of 
a laminar flow bench facility in the Engineering Depat tment was 
provided by Professor C. R. V i s i wanathan of the Electrical Sciences 
Division. 

When the modules were separated, it was observed that the 
headers of the fibre optic bundles were also coated with the same 
paint. Though the brass could be seen on the sharp edges, there 
was no flaking on large areas as there had been in the case of 
the detector housing. To remove these headers for cleaning and 
black oxide coating would have necessitated returning the instru- 
ment to GE for re-assembly and readjustment. The time delay would 
have forced cancellation of the flight program. 

Analysis of the data from the flights indicated the presence 
of loose flakes in the chopper compartment, particularly on the 
filters. The detector was subsequently removed and such pieces 
of dirt were observed as shown in Fig. 13- This represents a 
major but correctable defect, in the instrument prohibiting accurate 
analysis of the data from the flights. This effect has been 
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illustrated in the body of the report in Section 2.4.1. 

in re-assembl i ng Q the instrument, the cover of the mechanical 
module was rotated 45° to facilitate mounting in the instrument 
compartment of the scan platform. 

3.8 Additional C o mments and Conclusions . - Since proper 
determination of the calibration factors of the instrument is 
essemtial in data reduction and analysis, exhaustive discussion 
of this has been presented in the main report. (some arbitrariness 
has been exercised in deciding which particular items should 
be detailed in this appendix). To summarize, we will briefly 
reiterate major points concerning this as presented in both report 
and appendix: 

1) An elliptical distortion equivalent to 4% mixing of Q and 
U was observed durinq calibration factor determination prior to 

f 1 i ght. 

2) The character of this distortion indicates strongly 
that it is due to slight mi sor i entat i on of the polarizing prisms, 
but attempts to measure the orientation neither confirm nor deny 
this as a cause. 

3) Parameter mixing in rhe electronics has been definitely 
ruled out as a cause of the distortion. 

4) Analysis of the effect of the non-symmet r i c angular 
acceptance of the polarizing prisms shows that this cannot account 
for the elliptical distortion since such effects would introduce 
odd harmonic terms of the polarization angle in the T]-plane (Q/l, - 
U/l ). 


5) Until the cause of this distortion is pinned down there 
is some indefiniteness in the absolute value of the polarization 
determined by the instrument. 

6) This effect is, however, minimal in interpreting the 
flight data due to other problems, mainly the effects of the 
external windows and the loose dirt in the chopper compartment. 

7) It is felt that this distortion can eventually be 
accounted for and will present no significant problem. 

Going beyond this, we point out the major problems with the 
instrument which could be corrected: 

1) The remaining black paint must be removed from the fibre 
optics bundle headers and residual dirt removed from the chopper 
compartment. Also, the edge of the detector housing should be 
trimmed down to eliminate grazing by the chopper. 
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Fig. 16 -- Flakes of Black Paint on One of 
Interference Filters. This photograph was 
through the opening for the detector mount 
removing the detector from the instrument. 
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2) It would be desirable to replace the polarizing prisms 
with those with symmetric angular response, some additional 
expense mightibe incurred but there should be no reason why a 
more suitable cement cannot be found. T. Wise's observation that 
the optical path deviates from that expected indicates that these 
prisms are not as well fabricated as they should be. 

3) Some arrangement is desired to cope with the saturation 
problem (discussed in the body of the report). One such possi- 
bility is some real time indicator of saturation to enable the 
operator to lower the preamp gain manually. (The current optional 
gain is too low and in the cases encountered would have dropped 
the signals in the unsaturated colors to undesirably low levels.) 
Alternatively, it might be arranged to have different preamp gains 
for the first two and last two color channels. This might necessi- 
tate another sync signal generated by the filter wheel to elec- 
tronically switch gains. 

1) The 0.6 mv noise figure with barrels blocked as observed 
in flight represents achievement of the •§• millivolt error limit 
matching the drift limit of the processor electronics. 

2) The instrument and its electronics out-performed our ex- 
pectations. Comparison of the flight data with similar data from 
other pol ar irneters indicates the superiority of this instrument 
with regard to low scatter. 

3) The rate of data acquisition is virtually instantaneous 
compared to polarimeters typically used (though to be honest, there 
is technically no reason why rotating prism polarimeters could not 
be constructed with similar data acquisition rates). 

4) The instrument is sturdy and compact. 
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